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ABSTRACT 

Recently, it has been suggested that the metallicity aversion of long-duration gamma-ray bursts 
(LGRBs) is not intrinsic to their formation, but rather a consequence of the anti-correlation between 
star-formation and metallicity seen in the general galaxy population. To investigate this proposal, we 
compare the metallicity of the hosts of LGRBs, broad- lined Type Ic (Ic-bl) supernovae (SNe), and 
Type II SNe to each other and to the metallicity distribution of star-forming galaxies using the Sloan 
Digital Sky Survey (SDSS) to represent galaxies in the local universe and the Team Keck Redshift 
Survey (TKRS) for galaxies at intermediate redshifts. 

The differing metallicity distributions of LGRB hosts and the star formation in local galaxies forces 
us to conclude that the low-metallicity preference of LGRBs is not primarily driven by the anti- 
correlation between star-formation and metallicity, but rather must be overwhelmingly due to the 
astrophysics of the LGRBs themselves. Three quarters of our LGRB sample are found at metallicities 
below 12+log(0/H) j 8.6, while less than a tenth of local star-formation is at similarly low metallicities. 
However, our supernova samples are statistically consistent with the metallicity distribution of the 
general galaxy population. Additionally, we show that the star- format ion rate distribution of the 
LGRB and SNe host populations are consistent with the star-formation rate distribution of the SDSS 
galaxy sample. This provides further evidence that the low-metallicity distribution of LGRBs is not 
caused by the general properties of star-forming galaxies. Using the TKRS population of galaxies, we 
can exclude the possibility that the LGRB host metallicity aversion is caused by the decrease in galaxy 
metallicity with redshift, as this effect is clearly much smaller than the observed LGRB host metallicity 
bias over the redshift span of our sample. The presence of the strong metallicity difference between 
LGRBs and Type Ic-bl SNe largely eliminates the possibility that the observed LGRB metallicity bias 
is a byproduct of a difference in the initial mass functions of the galaxy populations. Rather, metallicity 
below half-solar must be a fundamental component of the evolutionary process that separates LGRBs 
from the vast majority of Type Ic-bl SNe and from the bulk of local star-formation. 



1. INTRODUCTION 

Shortly after long soft gamma-ray bursts (LGRBs) 
were identified as extragalactic events spanning cosmo- 
logical distances, it became apparent that they are fre- 



irregulars ( 


Fruchter et al. 


1999 


2006 


Le Floc'h et al. 


2003 


2002 . LGRBs have shown 


a strong preference for 



occurring m starforming galaxies (IFruchter et al. 


1999 


Christensen et al. 


2004 Le Floc'h et al. 20061), 


which 


often exhibit brig 


it emission lines (Bloom et al. 


1998 


Vreeswijk et al.|2001 Levesque et al.|2010a) indicative of 



have also frequently been associated with Type Ic super- 
nova (IStaiick et al.||2003| |Hjorth et aLl|2003| [Woosley fe 
Bloorh||2006 ) whose spectral lines show broadening from 
their high-velocity (~ 15,000 km s _1 ) ejecta. While 
a broad-lined Type Ic (Ic-bl) supernova (SN) has been 
found underlying the light of nearly every LGRB in 
which a deep spectroscopic search could be performed. 
The converse however, is most certainly not the case. 
Type Ic-bl SNe events without LGRBs are vastly more 
common than those with an accompanying LGRB. Thus 
LGRBs are now seen as a rare subtype of core-collapse 



supernovae (CCSNe) set along a progressive narrowing 
of stellar types which raise a number of astrophysical 
questions: What makes some CCSNe Type Ic's? What 
makes some Type Ic's relativistic broad-lined events, and 
what makes some of these relativistic broad-lined events 
also have a gamma-ray burst? The search for the phys- 
ical constraints that produce stars along this narrowing 
pyramid of phenomena is presently a primary focus of 
the study of massive stellar evolution. 



Several years ago, Fruchter et al. (2006) performed 
a detailed study of the LGRB host galaxy popula- 
tion, using the hosts of the Great Observatories Origins 
Deep Survey (GOODS) CCSNe sample as a comparative 
group. The results of this work showed a surprising dif- 
ference between the two populations. While half of the 
GOODs CCSNe occurred in grand design spirals, with 
the other half in irregulars, however only one out of 42 
LGRB host galaxies was in a grand design spiral. If one 
constrains the LGRB host population to a redshift of 1.2 
or less so as to match the redshift distribution of the 
GOODS SNe sample this drops to one out of 18, still 
a rather surprising result. The remainder of the LGRB 
host population are composed of generally faint, blue, ir- 
regular galaxies. This strong bias galaxy type has been 



2 



supported by later work using an even la rger comparison 
sample of CCSNe flSvensson et alpMol ). 
In additio n to the galaxy type preference, |Fruchter 



et al.| (120061) showed a strong preference for LGRB's oc- 
curring in the brightest, and hence likely the most star 
forming regions of their hosts. This suggests that LGRBs 
are formed from very massive progenitors (O and B stars) 
which often do not have time to travel far from their birth 
sites before exploding. However, CCSNe follow the blue 
light distribution of their host galaxies. This suggests 
that they do not generally require as massive progeni- 
tors as LGRBs. However, given that there is no evidence 
for a differe nce in the stellar IMF between blue irregulars 
and spirals ( Bastian et aLpOlO ), massive stellar progen- 
itors should be just as available p er unit star-formation 
in spirals as they are in irregulars 



(c.f . Tremonti et al.|2004 1 should be extended to a mass, 
metallicity, star-formation relationship. In this relation- 
ship the metallicity of a galaxy of a gi ven stellar mass 
is anti-correlated with its SFR. Indeed, Mannucci et all 



(20101 suggests that this single relationship may be able 
to extend to quite high rcdshift, with the generally lower- 
mctallicities of galaxies at higher redshift corresponding 
w ell to their higher SFR. 



Kocevski fc West| ( |2011[ ) and |Mannucci et al| ( |2011[ ) in- 
dependently argue that the apparent LGRB preference 
for low-metallicity hosts is thus due to a more funda- 
mental mass, metallicity, & SFR relation. LGRB hosts 
are low-metallicity because they are effectively selected 
based on the basis of star-formation. However this does 



not ex plain why, as already noted by |Modjaz et al. 



therefore concluded that LGRB formation likely requires 
a low metallicity progenitor. The bias towards irregular 
galaxies is then a result of their generally low metallicity 
as expected by the mass-metallicity relation. 

A s imilar conclusion was reached by |Stanek et al.| 
(2007), who showed that the very nearest LGRB hosts 
all have low metallicity when compared to similar mag- 
nitude galaxies in the Sloan Digital Sky Sur vey (SDSS) 
sample. Furthermore, Kewley et al. (2007) found the 
LGRB host sample to be comparable to extremely metal- 
poor galaxies in luminosity-metallicity relation, star- 
for mation rate (SFR), and i ntern al extinction. Addition- 
ally, Wolf fc Podsiadlowski| ( |2007| ) formally demonstrated 
that a metallicity dependence could indeed explain the 
difference in morphological type between the hosts of 
LGRBs and the more general CCSNe population. 

The suggestion that LGRBs form preferentially in low- 
metallicity hosts, however, is not unifo rmly accepted- 
Using a wide redshift sample of LGRBs, Savaglio et al. 
(2009) argued that LGRB hosts lie on the sam e mass- 
metal l icity r elation as regular galaxies. Indeed, |Berger| 
|et al| ( |2007[) have used this claim to argue that because 
the host of LG RB 020127 is unusuall y bright it must also 
be metal rich. jPeeples et al.| ( |2009[ ) however, show the 
existence of low metallicity outliers on the luminosity- 
metallicity relation, and argue against the assignment of 
metallicities (to individual galaxies) based only on their 
luminosities. In particular they highlight the morpholog- 
ically simila rities of their bright o utliers to the brighter 
hosts in the Fruchter et al. ( 2006 ) s ample. 



Fruchter et al.|(|2006|) (120081), the Type Ic-bl SNe without associated LGRBs 



In a comparison of CCSJNe hosts, Modjaz et al. (20081 
showed that LGRBs are observed to occur m host galax- 
ies with much lower metallicities than either the hosts 
of Type Ic-bl SNe or the bulk of the star-forming galax- 
ies in SDSS. Given the dramatic difference in metallicity 
between the Type Ic-bl and LGRB samples, which per- 
sisted even when hos t s of s imilar luminosity were com- 
( |2008[ ) concluded that this is due 



pared, Modjaz et al. 



to a metallicity avoidance among the LGRBs. As Type 
Ic-bl SNe are frequently observed to be coincident with 
many LGRB events this seemed to suggest a metallicity 
dependent step in either the formation of the requisite 
gamma-ray jet or in its ability to escape their progen- 
itor during a highly relativistic core coll apse SNe of an 



already ou ter layer stripped massive star ( Langer & Nor- 
~|[2006 ). 



man 



Recently however, Mannucci et al 
gested that the well know mass metal 



do not show a preference for low-metallicity hosts. Nor 
does it explain why CCSNe in the GOODs sample should 



Fruchter et al. 


2006 


Svensson et al. 


2010 



One would 



(2010) has sug- 



icity relationship 



star-forming galaxies, and thus low-metallicity. 

In order to more clearly address this issue we compare 
the metallicity distribution of the hosts of LGRBs with 
that of the hosts of several similar indicators of star- 
formation: LGRBs, Type Ic-bl, and Type II SNe. Fi- 
nally, we compare all of these hosts galaxies with more 
general catalogs of star-forming galaxies, compiled from 
the SDSS and Team Keck Redshift Survey (TKRS). This 
broad-range of star-formation indicators gives a good and 
varied base from which to determine whether LGRBs arc 
indeed biased toward low metallicity even in excess of 
that expected due to the high SFRs of their hosts. 

In the next sections we describe in detail the samples 
used for these comparisons. We then present the meth- 
ods employed to place all of the samples on a consistent 
metallicity scale. Finally we present the direct compari- 
son of the metallicities of these populations. 

2. THE SAMPLES 

The fundamental question of this paper is whether 
LGRBs track star-formation, or whether their choice of 
host is biased by metallicity. In order to determine this, 
we will want to compare the locations of LGRBs with 
other tracers of star-formation. 

The most obvious tracer of star-formation to compare 
with LGRBs are star-forming galaxies them selves. We 
use the Slo an Digital Sky Survey (SDSS - Abazajian 
et al| 2009 ) to obtain a sample of star- forming galax- 
les at lo w-redshift and the T eam Keck Redshift Survey 
(TKRS - |Wirth et aL]|2004[ ) to allow us to extend our 
galaxy star-forming sam ple o ut t o z ~ 1. These samples 
are described in Sections 12. 1| and |2.2[ respectively. 

Type II SNe also provide a very direct estimate of star- 
formation. While one might worry that some Type II 
SNe are lost to extinction, the same is true of the Ha 
emission which is used to estimate the star- format ion of 
galaxies and of (to a lesser extent) LGRBs themselves (a 
radio or X-ray position is often sufficient to identify the 
host of a LGRB). 

Many LGRBs are associated with Type Ic-bl SNe. In- 
deed in almost every case where one estimates that one 
should have been able to discern the light (or in quite 
a few cases the spectrum) of a Type Ic underlying the 
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afterglow, were it there , evidence for a Type Ic has been 
found (see Cano 2012 for a detailed discussion of this). 
However, where spectral observations have been done, 
they reveal Type Ic with broad-lines due to a very en- 
ergetic explosion imparting large velocities to the ejecta. 
Therefore, we also compare the locations of LGRBs with 
Type II and Type Ic-bl SNe, which are not known to be 
associated with LGRBs. 

Finally, we present our sample of LGRBs we have as- 
sembled using both our own spectra and spectra available 
in the literature. 

Below we describe all our samples in detail. 

2.1. Star-forming Galaxies from the SDSS 

The SDSS is a combined photometric and spectro- 
scopic survey of approximately 14,555 square degrees 
with spectroscopy over 930,000 galaxies. The SDSS gives 
accurate astrometry, photometry and line fluxes. The 
survey has already been used for numerous stu dies, and 
in particular the work of Tremonti et al. (2004) examin- 
ing the distribution of metallicities ot the galaxies. Here 
we will use this catalog as the basis of our comparison, 
but we rederive metallicities using the R23 method (as 
for the LGRB sample). Our metallic ity d erivation will 
be discussed in more detail in Section [3TTI 

The SDSS provides a pre-selected sample of star- 
forming g alaxies using the SUBCLASS starforming or 
starburst (Sanchez Almeida et al. 2010). These galaxies 
have lines indicative of star-formation, but do not appear 
to be dominated by an Active Galactic Nucleus (AGN). 
However, we further restrict these galaxies based on line 
strengths in order to assure we can obtain reasonably 
accurate metallicities. 

As described in Section |3.1.1| the R23 diagnostic re- 
quires the 3727 A [O II] doublet, the 4861 A H/3 line 
and the 4959 A & 5007 A [O III] lines. The 3727 A 
[O II] doublet, commonly used by the R23 method, is 
actually a combination of the 3726 A and 3729 A lines, 
which are frequently unresolved in long slit spectroscopy. 
Here these lines are summed and their errors combined 
statistically to form the 3727 A line before applying the 
signal to noise cut in order to better replicate the long 
slit methodology. The [O III] lines are in a 3:1 flux ra- 
tio; therefore we do an optimal combination to generate 
improved line fluxes and respective error before apply- 
ing the signal to noise cut. In both cases independent 
errors between lines are assumed. The 4861 A H/3 and 
6563 A Ha lines are processed with a simple S/N > 8 
cut. No similar constraints are placed on the 6583 A [N 
II] line to avoid introducing a metallicity bias (i.e. under 
detection of low metallically systems due to discarding 
galaxies with low [N II] S/N). Since the 6583 A [N II] 
line is in such close proximity to Ha, the Ha reductions 
are sufficient to remove bad cases on [N II] . Also with the 
exception of the [N II] all specified lines must be nonzero. 
It should be noted that the requirement for detection 
of the 3727 A line and the limits of the SDSS spectral 
wavelength coverage limit the comparison population to 
a redshift of 0.0209 or higher. 

The SDSS catalog we have actually employed is the 
Max-Planck-Institut fur Astrophysik - John Hopkins 
University (MP A-JHU) emission line analysis for SDSS 
Data Release 7 ( |Brmchmann et aL]|2004 Tremonti et al. 



2004^ In order to fully reproduce our sample again 
the user should specify SPECTROTYPE galaxy, SUB- 
CLASS starforming or starburst, and also require non- 
zero ugriz CMODEL values, which we required for sub- 
sequent photometry. The resulting sample consists of 
approximately 137 thousand galaxies. 

From our full SDSS sample, we also create two sub- 
selected samples. The first employs a redshift cut to 
create a volume limited sample. This is done to minimize 
the incompleteness of faint galaxies in the sample. The 
second sample is a weighted sample of the first, where 
the weighting chooses galaxies as SNe and LGRBs may 
do: by their SFR. 



2.1.1. 



redshift of z < 0.04 SDSS general star- 
galaxies 



-forming 



The SDSS is a magnitude limited sample. This means 
that large galaxies (generally at significant redshifts) are 
overrepresented compared to what one finds in a fixed 
volume. This has the effect of skewing the metallic- 
ity distribution, due to the previously mentioned mass- 
metallicity or nearly equivalently, luminosity-metallicty 
relationship of galaxies. To obtain an unbiased sample 
of metallicities a volume limited sample is preferable. 

We are forced to consider only galaxies with a redshift 
greater then z ~ 0.02 by our requirement of the 3727 A 
[O II] doublet for the R23 method. At lower redshifts, 
the doublet falls to the blue of the SDSS spectral range. 
Our ability to reach high redshifts is limited by the mag- 
nitude limit of the SDSS survey. The SDSS spectroscopic 
survey is only complete to m 9 = 18, and becomes very 
sparse at fainter magnitudes (some fibers, but not many, 
were devoted to objects fainter than m g = 18). The ef- 
fect of this magnitude cut on the available comparison 
SDSS sample as a function of absolute magnitude can 
be seen in Figure [T] Similarly, in Figure 12] we show the 
distribution of star-formation in our SDSS sample as a 
function of metallicity for several redshift cuts. A min- 
imal redshift range extending 0.02 < z < 0.03 would 
give completeness down to a luminosity of ~ -17.5, 
however our ability to constrain the metallicity distribu- 
tion of galaxies brighter than Mg ~ -21 is limited by a 
small number statics. Extending the redshift range to z 
< 0.04 gives a sufficient sample size on the bright end 
and completeness down to a luminosity of Mb ~ -18. 
Further extending the redshift range to z < 0.05 gives 
no significant improvement on the bright end and would 
reduce the completeness down to a luminosity of Mb ~ 
-18.5 and is thus not advantageous. Our resulting red- 
shift range is identical to that chosen by |Prieto et al.| 
( |2008[ ) in their supernova host survey for similar reasons. 
This leaves approximately 21 thousand galaxies within 
the 0.02 < z < 0.04 redshift cut. 

If we fit a power-law to the star-formation function 
shown in Figure [T] from Mg ~ -19 to the completeness 
luminosity limit or redshift range of Mb ~ -18, we find 
a slope of log(SFR) cx -0.3 x Mb- Our star-formation 
numbers are derived from Ha luminosities, but if one as- 
sumes that star-formation is proportional to galaxy lu- 
minosity, this corresponds to a Schechter function slope 
of a = —1.3, or equal to that derived by |Blanton et al.| 



1 See http://www.mpa-garching.mpg.de/SDSS/ for the data 
products, their descriptions, and a more detailed citations list 
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(|2005[). This ag rees well with that determined by Blan- 



ton et al. (2005). It implies that approximately 21% of 
star-formation lies in galaxies with an Mb dimmer then 



-18. However, Blanton et al. (20051 argues that the ac- 
tual slope at these tainter magnitudes might actually be 
as steep as a = —1.5, in which case approximately 31% 
of star-formation lies in Mb > -18 galaxies. The num- 
bers are in contrast with our incomplete sample with only 
about 6% of star- format ion occurring in galaxies dimmer 
than an absolute magnitude of -18. While in many cases 
we limit our discussion to galaxies brighter than our - 
18 effective limit, but where we do not, these fractions 
should be kept in mind. Imposing an Mb brighter then 
-18 magnitude cut on the volume limited SDSS sample 
yields approximately 15 thousand galaxies. 

2.1.2. Synthetic star-formation weighted galaxy sample 

The SDSS selects galaxies by number. All galaxies 
large or small, if detected, appear as a point on the 
SPSS plot. However, SNe and LGRBs (if, for example, 
|Mannucci et al.J|2011| are correct) , choose their galaxies 
with a probability proportional to the star- format ion of 
the host. It is therefore useful in this paper for illus- 
trative purposes to create a sample of galaxies weighted 
according to star-formation. To do this, we take our 
SDSS volume-limited sample, and select galaxies from it 
with a probability of selection proportion to the SFR of 
the galaxy. In practice this means that in our random 
weighted sample, a given galaxy from the original sample 
may not appear, or it may appear more than once. 

2.2. z < 1 Star-Forming Galaxy Population Via TKRS 

Our desire to obtain a volume limited galaxy sam- 
ple with a good representation of faint galaxies has lim- 
ited our comparison SDSS sample to very low redshift 
(z < 0.04). However, in order to get a reasonable sample 
of LGRBs, we include all LGRBs with measured emission 
line metallicities. While this sample has a median red- 
shift of 0.3 (the effective extent of the full SDSS sample), 
it extends up to a redshift of 0.8. Therefore if we are to 
also take into account the evolution of galaxy properties 
with redshift in our comparison, the SDSS sample is an 
inadequate sole comparison. 

While no available sample perfectly meets our require- 
ments, the TKRS met allicity sample provided in Kobul- 



nicky fc Kewley| ( 2004 ) is excellent in many respects. The 
sample has an effective magnitude (or more appropriately 
line-strength limit) similar to that of the LGRB sample 
that we have obtained. It extends out to z ~ 1, and thus 
agrees well with the reds hift range of our LGRB sam- 
ple, and is already on the Kobulnicky & Kewley (20041 



metallicity scale, which we use. 

The primary drawback of the TKRS sample is that 
the [N II] /Ha observations necessa ry to r esolve the i?23 
metallicity degeneracy (see Section 3.1.1) are not avail- 
able in all cases - primarily on galaxies at z 0.5 where 
the [N II] and Ha lines move into the infrared. In these 
cases, where the branch of the R23 metallicity diagnos- 
tics is uncertain, the TKRS sample assumes the upper 
branch. This is a clear bias. However, at redshifts z ^ 0.5 
only about 10% are in the lower branch. While higher 
redshift galaxies of similar luminosity tend to move to 
lower metallicities, this magnitude limited sample also 



moves to higher luminosities at higher redshifts. Thus it 
is reasonable for us to assume that at redshifts above 0.5 
perhaps 10% of the sample, all in the lower branch, are 
missed. 

2.3. Supernovae Host Galaxy Populations 

In addition to the LGRB hosts and the low redshift 
general star-forming galaxies, we also include supernovae 
(SNe) host galaxies in our analyzes and resulting com- 
parisons. SNe are direct tracers of the deaths massive 
stars - and due to the short lifetimes of massive stars, 
their deaths trace star-formation. Broad-lined Type Ic's 
(Ic-bl's) are often found under the decaying afterglow 
of an LGRB. Thus studying the metallicity distribution 
of these explosions provides a direct test of whether the 
metallicity bias that may be found in LGRBs is only 
seen in LGRBs, or whether it is also shared by a larger 
class of particularly energetic exploding stars. However, 
in order to provide an unbiased comparison sample we 
supplement the Type Ic-bl SNe with the more readily 
available Type II events. The majority of SNe are of 
Type II and these are perhaps the archetypical SNe. Re- 
cently it has been shown that not only do these SNe 
track blue light between galaxies; their format ion loca- 
tions also t r ack the blue light within galaxies (Fruchter 



et al. 2006 Kelly & Kirshner 2012). They may there 



tore be the most-unbiased tracers of star-formation be- 
sides the massive stars themselves. For maximum consis- 
tently in both data and subsequent analysis, we use only 
Type II whose hosts are in our volume limited SDSS 
star-forming galaxy population. This is a luxury that 
the relative frequency of Type II SNe allows. Comparing 
between the SNe types allows us to be sensitive to poten- 
tial biases introduced by the broader redshift range for 
the broad-lined Ic's or the different data sourcing for the 
Type Ic-bl population as well as any intrinsic differences 
between the SNe types. 

2.3.1. Broad-lined Type Ic SNe hosts 

For our broad-line d Type Ic (Ic-bl) sam ple we build 
upon the sample of Modjaz et al.| ( |2008[ ) and expand 
it with a few additio nal objects that are now available 
Sand ers et al.||2012|). As addressed in Modjaz et al. 



2008p biases are introduced when SNe are found via 



targeted surveys. In these surveys, such as the Lick Ob- 
servatory Supernovae Search (LOSS) with the Katzman 
Automatic Imaging Telescope (KAIT), a galaxy popula- 
tion pool is selected in advance then repeatedly searched 
for the occurrence of SNe. This obviously has potential to 
bias the SNe population with the selection criteria used 
in compiling the galaxy search pool. Since this has the 
potential to affect the metallicity luminosity relation of 
those SNe we distinguish between SNe found in targeted 
and untargeted surveys whenever displaying the objects 
in plots. 



Also addressed in Modjaz et al. (2008) are differences 
between the metallicity at the SNe location with in the 
host and the metallicity at the host center. Since in 
most cases we lack mea surements of both, w e follow the 
general methodology of Modjaz et al. (2008) employing 
a metallicity gradient to estimate the metallicity at the 
SNe location when it was not specifically measured. (See 
Appendix [A] for details on how this fit was determined 
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Figure 1. Histogram plots of the number (left) and star- formation (right) distribution of the SDSS population at various redshift ranges. 
Note that there is a lack of any significant star-formation contribution in galaxies dimmer then -17 Ms and this contribution drops off 
around -18 Ms for the redshift cut populations shown. For the 0.02 < z < 0.04 redshift cut adopted we fit star-formation in the -18 > Mb 
> -19 range and extrapolate to esti mate the Mb > -18 star-formation using both the slope as fitted and the slope as estimated from the 
c?2 ~ —1.5 value estimated inlBlanton et al.l 120051. 
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Figure 2. Fractional cumulative distribution plots of star-formation vs. metallicity for various redshift cuts in the SDSS general star 
forming galaxy population. From this it is apparent that the volume limited (0.02 < z < 0.04) SDSS population is complete only for objects 
brighter then -18 M^. 



and applied). In the case of SN 2002ap, occurring in M74 
/ NGC 628, we do not need to estimate the metallicity 
gradient of the ho st. Instead we adopt spectra from |Fer-"1 
guson et al. ( 1998 ) using the closest spectra from the core 
to calculate the central galaxy metallicity and estimate 
the SNe site metallicity from spectra taken at a similar 
radius from the galaxy core. 

2.3.2. Type II SNe hosts 

A large s ample of Type II SN e hosts are available from 
the work of |Prieto et al. (2008). Limiting the sample to 



Type II SNe within the redshift range outlined in Section 
1 2. 1.1 1 and hosts present in the SDSS general star- forming 
galaxies sample yields 16 objects. However 7 of these 
SNe were discovered via targeted galaxy surveys and thus 
may potentially be biased to the brighter host galaxies 
(as seen in the targeted Type Ic-bl population). In an 
attempt to track this bias we also employ a population 
of non-targeted Type II SNe for comparison. 
To expa nd our untargeted Type II sample, we use 



the SDSS ( [Frieman et al.||2008| |Sako et al.||2008D and 



Supernova Factory (SJN Factory - Aldering et al. ' 2002 1 



G 



SNe searches selecting only Type II events with hosts 
in our volume limited SDSS star-forming galaxy popula- 
tion. Unfortunately, this yields only 3 and 7 objects from 
the SDSS and SNFa ctory surveys respe ctively (of which 4 
were already in the Prieto et al. |2008| SNe sample). Fur- 
thermore, there is evidence that these surveys are biased 
towards preferentially detecting events on low surfaces 
brightness backgrounds as described in Appendix |B| 

2.4. LGRB host galaxies 

To provide as comprehensive an LGRB host galaxy 
population as possible we co mbine our own s pectroscopic 
observations of LGRB Hosts (Graham 2012 in prep) with 
those from the literature. Our combined sample contains 
all hosts known to us which have emission line spec- 
troscopy capable of giving a metallicity value using the 
R23 method. When using observation from the literature, 
we compute our own metallicity values via the published 
line strengths rather then adopt published metallicities 
so as to obtain better consistency between samples. The 
line strengths used as well as the derived quantities such 
as calculated metallicity and SFR are given found in Ta- 
bles [TJ [2] & [3] presented in Section |3.3| 

The metallicity value of the location of a LGRB or SN 
will generally not be exactly that of the center of the 
host. However, LGRBs frequently occur at or near the 
brightest regions of their generally irregular host galaxies 



(Fruchter et al. 2006 Svensson et al. 20101. Therefore, 



when a spatially resolved spectrum ot a LGRB host is not 
available, it is reasonable to use the average host value 
as the value for both the galaxy center and the site of 
the explosion. 

There are however three LGRBs in our sample where 
the LGRB site region and host galaxy center are well 
separated at ground-based resolution: LGRBs 980425, 
020819B, 060505. These deserve special treatment. 

1. LGRB 980425/SN 1998bw is the nearest known 
LGRB. It lies on a dwarf irre gular, but because of 



its low redshift (z ~ 0.008 - Galama et al. 19981 



the host is very well resolved from the ground. 
Ther e are several published spectra of the LGRB 



site (pollerman et al. 2005 Hammer et al. 2006 



Christensen et al. 2008 1 from which we determine 



metallicity value s of lo g(Q/H)+12 = 8 .48 p olled 
man et ai.||2005[), 8.64 |Hammer et al ||2006| , and 



8.23 ( CEristensen et al.p008|). We adopt the me- 



dian value, 8.48 from Sbllerman et al. 

We note, however. 



p05] 

that a 



for 
dis- 
II] di- 



the site metallicity. 

agreement betwe en the R23 and [Nil] / [O 
agnostics for the Sollerman et al. ( 2005 1 spectrum 
suggests that the errors on that spectrum may be 
large. 



The IFU data of Christensen et al. (2008) allows us 
to make a rough estimate of the central metallicity 
of the host. The fibers in the SDSS, subtend on a 
scale of ~ 1.2 to 2.3 kpc in our volume limited sam- 
ple. Sampling the IFU elements within radii 7.4" 
and 14.2" of the core, we compute metallicities of 
8.59 and 8.55 respectively. If we take our LGRB 
site metallicity value and apply the metallicity gra- 
dient method described in Appendix |A"| to estimate 
the central metallicity we also find a value of 8.55 



which we adopt as our central metallicity value for 
LGRB 980425. 



2. The host of LGRB 020819B is a large spiral ( |Jakob- 
sson et al.||2005 Graham 2012) with a diameter of 
about 7" on the sky. The LGRB was dark, but a 
radio position with an error circle of 1" diameter 



( |Jakobsson et al.||2005 ) places it on a small "blob" 
at the edge of the spiral (this could either be an 
outlying star-forming region of the spiral or pos- 
sibly a satellite). The metallicity of the center of 



the ho st has been well determined ( Levesque et al 



2010b). However, a problem with the blue chip ot 



the spectrograph on the same observing run, meant 
that a full R23 metallicity could not be obtained for 
the "blob" under the location of the burst. How- 
ever, the [N II] /Ha diagnostic can be us ed to give 



a rough estim ate of the local metallicity ( Levesque 
et al.|2010b[). W e will discuss this object 
Graham"! 2012) in prep 



urther m 



3. The host LGRB 060505 (z ~ 089) is a well re 



solved spiral ( |Th5ne et al |2008| . The LGRB oc- 
curred right on top of an Hli region. No supernova 
was seen and strong limits have been plac ed on 



the magnit ude of the underlying supernova (Ofek 
et al. 1120071) . There has been some discussion as to 



wheth er "this burst was actually a long or a short 



burst (Ofek et al 
ever, due to the 
region 



2007 



2008). How- 



Thone et al. 
ocation of the burst on an HII 
and the bursts ~ 4s duration, we have in- 



cluded it in our long burst sample. |Th5ne et al.| 
( 2008 ) gives us metallicities for both the center of 
the host (Bm in their table) and burst location. 



3. METHODS 



DETERMINATION OF ANALYZED 
PROPERTIES 



The majority of analyses used within this paper are 
based on values of four physical galaxy properties: metal- 
licity (both event local and galaxy central), total SFR, 
absolute B band magnitude (Mb), and redshift. Here 
we describe how these values (except for redshift) are 
determined. 

3.1. Metallicity 

Due to differences between metallicity diagnostics and 
their various calibrations, a true comparison of metal- 
licity requires using a common scale and if possible a 
con sistent methodology an d diag nostic. Here we adopt 
the Kobulnicky & Kewley (2004) scale exclusively and 
the R23 diagnostic to the extent possible. To ensure 
consistency of method and scale we calculate metallic- 
ity va lues ourselves from published fluxes and our own 
work ( Graham|2012| in prep). Metallicity calculations are 
preform ed with an improved versi on of the IDL cod e out 



lined in Kcwlcy & Dopita (2002) (updated to the Kob 



ulmcky & Kewley||2004| scaleJT 



dr conversion to other 
scales and discussion of associated issues we refer the 



reader to Kewley & Ellison (2008) and the metallicity 
diagnostic transformati ons given there in. 



Based on the |Allende Prieto et al.| ( |200l[) 6300 A [O 
I] line measurements of the Sun, solar me tallicity is esti- 
mated to be log(Q/ H)+12 = 8.69 ± 0.05 (|Allen d e Prieto j 
et al.||2001[ ) in the |Kobulnicky fc Kewleyj ( |2004[ ) scaleT^ 
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It should be noted that converting solar metallicities to 
emission line H II region abundances is exce edingly diffi- 
cult. Thus the |Kobulnicky fc Kewley||2004| scale is more 
accurate internally than to any absolute reference to solar 
value and such comparisons should be avoided or limited 
to broad generalizations. 

3.1.1. R23 Diagnostic 

The R23 method is a commonly used metallicity di- 
agnostic based on the electron temperature sensitivity of 
the oxygen spectral lines, achieved using the ratio of oxy- 
gen line strengths to a hydrogen spectral feature, which 
is independent of metallicity . Doubly ionized oxygen, [O 
III] , has strong lines at 4959 and 5007 A and [O II] has a 
particularly strong, typically unresolved, doublet at 3727 
A. The metallicity independent 4861 A H/3 line is conve- 
niently placed between these sets of oxygen lines. This 
close placement minimizes errors due to misestimates of 
internal extinction and possibly even systematic effects 
of equipment. The ratio of the sum of the fluxes of these 
oxygen lines divided by the H/3 flux gives the equation for 
R23 used in the classical application of this diagnostic. 



R-2 



I: 



3727 



'4959 



5007 



H/3 



(1) 



First pro posed by Bernard Pagel in 1979 (Pagel et al. 
1979 1980), R23 has become the primary metallicity di- 
agnostic for galaxies where the faint [O III] 4363 A line is 
not measurable, and thus has been used extensively for 
work at moderately high redshifts, where the majority of 
our LGRB sample is found. 

In order to correlate the flux of a line belonging to an 
individual atomic ionization level with the total abun- 
dance of that element it is necessary to know what frac- 
tion of the element is ionized to the level in question. 
Here, in the case of oxygen, this is achieved by measur- 
ing the flux ratio between the [O II] and [O III] lines. 
This gives the relative population in the O II and O III 
ionization states, and allows one to fit the metallicity for 
that specific ionization state ratio. 

Thus in the classical application of this diagnostic, the 
R23 value would be calculated from the measured line 
ratios and then compared along an [O III] to [O II] line 
ratio contour. This classical application, however, treats 
ionization as a parameter independent of metallicity and 
ignores the effect metallicity has on t he ionization state 
of the lines. Kewley & Dopita ( 2002[ ) solve this issue by 
iterative fitting of the metallicity and ionization param- 
eter. 

While we have so far mentioned only oxygen line fluxes, 
in some cases observational difficulties may mean one has 
only equivalent line widths for a g iven spectrum. How- 
ever, [Kobulnick^~^^hilIips ( 2003[ ) has established that 
the R23 method can be directly used on equivalent width 
values (instead of reddening corrected flux values) and is 
found to be more accurate than flux ratios when redden- 
ing information, and in particular the Balmcr decrement, 
is not available. 

The R23 method diagnostics suffer a degeneracy due to 
different effects being dominant at different regimes. In 
the low metallicity regime the effects of the metals on the 
electron temperature of the system can be ignored due 
to their low relative abundance. Thus the more metals 



in the electron gas, the more collisional excitations and 
more resultant flux in the metal lines. As the metallicity 
rises however, emission from infrared fine-structure lines 
becomes significant and serves as a cooling mechanism, 
lowering the electron temperature, the electron velocity, 
the number of collisional excitations, and thus the metal 
line flux. 

This temperature dependence causes two metallicity 
values (one high, one low) to generate the same R23 line 
ratio. Unless one of the degenerate values can be obvi- 
ously excluded (or the two values are within the error 
range of the R23 calibration) new empirical data is the 
only accurate way to break the degeneracy. This is usu- 
ally accomplished, as we do here, by application of the 
[N II] /Ha diagnostic. 

3.1.2. [N II] /Ha Diagnostic 

The [N II] /Ha line ratio provid es a crude metallicity 
indicator (Kewley & Dopita 2002). However, due to its 
strong dependence on the ionization parameter it gives 
only a gross estimate of abundance unless the ionization 
is known. Since determining the ionization parameter 
requires measurement of lines which themselves consti- 
tute better metallicity diagnostics, the application of an 
ionization parameter correction to this diagnostic is of 
limited utility. Also the diagnostic is easily distorted 
by contamination from shock excitation or hard ionizing 
radiation from AGN. Thus its primary application is se- 
lecting between the degenerate upper and lower branch 
values provided by the R23 diagnostic. 

However use of the Ha and 6583 A [N II] lines can be 
advantageous especially when a low redshift coupled with 
instrumentation limitations preclude measurement of the 
3727 A [O II] line. Additionally the small separation 
between the 6584 A [N II] line and Ha lines renders this 
method immune to the effects of reddening and enables 
it to be applied without any flux calibration required. 

While the degeneracy of the R23 diagnostic can be bro- 
ken, as we do in this work, by application of the [N II] /Ha 
diagnostic, this has the obvious disadvantage of requiring 
additional measurements considerably outside the R23 
wavelength range This in turn often necessitates a sepa- 
rate observation and for higher redshift objects use of a 
separate near-infrared, instrument. Nonetheless, the [N 
II] /Ha diagnostic is one of the most reliable methods for 
breaking the R23 degeneracy. 

3.2. Star- Formation Rates 

We estimate S FRs from galaxy Ha fluxes using the 
Kennicutt ( 1998[ ) Ha SFR diagnostic. When an Ha flux 
is unavailable (for instance a LGRB host above z = 0.6 
without calibrated infrared observations), we estimate 
the Ha flux using the H/3 flux via the Balmier decre- 
ment). Unlike metallicity, estimation of SFRs requires 
absolute (as apposed to relative) line flux measurements 
over the entire galaxy. This requirement significantly 
complicates our analysis, as we are typically dealing with 
slit or fiber spectroscopy which does not cover the entire 
galaxy. Therefore, in addition to a standard absolute 
calibration of the spectroscopy, we must correct for slit 
or fiber losses. We discuss these corrections below. 

3.2.1. Slit loss correction for LGRB host galaxies 
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To obtain slit loss estimates for the LGRB popula- 
tion in cases where the host is not well-resolved from 
the ground, we convolve HST images of the hosts to 
ground based seeing and then determine the fraction of 
blue light that would land in a slit sized box. When pos- 
sible these slit loss fractions were computed for the size 
of slit used (typically 1" and assumed to be 1" when not 
stated explicitly) after convolving the image to match 
the described seeing (or typical seeing for the telescope 
used when not stated explicitly) in two orientations (one 
maximally the other minimally favorable). 

Typical slit loss estimates were between 25 to 40 % 
with a typical difference between the orientations of un- 
der 5%. Given the small size of the orientation depen- 
dance the average between the orientations is used for the 
correction. Due to the stability of the flux loss, in cases 
where this correction can not be computed (i.e. sufficient 
images are not available) the average value is used. This 
approach does have the obvious limitation of assuming 
that the measured line strength is representative of the 
object as a whole (central galaxy spectroscopy is used 
when possible). Cases without suitable spectroscopic es- 
timate of star-formation are dropped from analyses that 
require the actual host SFR; otherwise the entire sample 
is used 

3.2.2. fiber loss correction for SDSS galaxy derived 
populations 

Just as with slit spectroscopy, it is necessary to correct 
for the difference in total flux of the SDSS objects with 
the flux observed propagating down the fibers. Fiber loss 
estimates are computed using the flux ratio of the fiber 
plug magnitudes to the estimate of the total galaxy mag- 
nitude, which is reported as CMODEL magnitudes in 
the MPA -JHU emission line analysis (described in Sec- 
Galaxies without CMODEL magnitudes are 



2.1) 



tion 

discarded. This methodology does have potential issues 
with biasing the galaxy SFR to the center of the galaxy 
intensity. However while this source of error can not 
be effectively eliminated, is also a slight issue with the 
LGRB population as well, and is countered as best as 
possible by using the galaxy central properties for com- 
parison. Since the Type II SNe hosts samples as well as 
the synthetic star- format ion waited random selection are 
both subsets of the SDSS spectroscopic sample, they are 
corrected as well as part of this procedure. 

3.2.3. SDSS spectroscopy substitution for Type Ic-bl SNe 

hosts 

The Type Ic-bl SNe hosts population is highly inho- 
mogeneous. The orientations and sizes of the slits are 
often not well documented, and due to the the proximity 
of these hosts, the slits subtend very small fractions of 
the galaxies. Due to the absence of accurate slit loss es- 
timates on the Type Ic-bl host spectroscopy and the low 
redshift / large angular size / low slit coverage / high slit 
loss of this population it is not directly us able in analyses 
that require the actual host SFR (Sections 4.3.2 & Figure 
[6|. However rather then omit the Type Ic-bl SNe hosts 
completely we search the SDSS sample for spectroscopy 
on their host galaxies and use that to calculate SFRs as 
in the previous section. As this spectroscopy is used only 
for SFRs (and not metallicities) coverage of the 3727 A 
[O II] line is not required thus allowing application of this 



methodology without the z > 0.0209 redshift limitation. 
This yields 7 Type Ic-bl SNe hosts with usable SFRs. 

3.2.4. M B 

For the SDSS galaxy derived populations absolute 
B band magnitudes are determined from the redshift 
and CMODEL ugriz magnitudes provided in the MPA- 
JHU-SDSS data products via the SDSS2BESSEL L rou- 



tine of the KC ORRECT IDL software package (Blan- 
ton & Roweisl 120071). The provided CMODEL rnagm- 



tudes are integrated galaxy apparent magnitude esti- 
mates for the spectroscopic targets based on the SDSS 
photometric data. The SDSS2BESSELL software pack- 
age converts the ugriz magnitudes into the UBVRI mag- 
nitude system (and requires input values in all 5 SDSS 
bands). Conversion from apparent to absolute magni- 
tude using the SDSS redshift value is preformed using 
the LF_DISTMOD IDL software package. 

Mb values for the LGRB host galaxies and Type Ic-bl 
SNe are taken from the literature (sometimes as apparent 
magnitudes and converted to absolute via redshift, in the 
same manner via LF_DISTMOD, as needed) . Mg values 
for the TKRS sa mple are taken directly from |Kobulnicky| 



fc Kewley| ( |2004| (along with all other TKRS data used). 

We choose to employ the luminosity metallicity rela- 
tion in place of the mass metallicity relation primarily 
due the availability of B band luminosity data across 
our populations. Many of our objects, both LGRBs and 
Type Ic-bl SNe hosts, do not have the necessary pho- 
tometric coverage to allow for a robust determination 
of mass values. Attempting to extrapolate photometric 
coverage from spectroscopic continuums is difficult given 
the faint continuums present in most of our emission line 
dominated spectra, especially of the LGRB host galax- 
ies. The Type Ic-bl SNe also have the difficulty of being 
close enough that long slit spectroscopy gives only lim- 
ited spacial coverage of the host (which in many cases 
is not at the galaxy center) thus the colors would likely 
not be indicative of an host galaxy as a whole. Luminos- 
ity also has the added benefit of avoiding coverage issues 
and biases at different redshifts (i.e. ugriz coverage with 
SDSS gives very different rest frame coverage at z ~ 
vs. z ~ 0.8 while the B band luminosity can be easily ex- 
trapolated througho ut our entire redshift range). Since 
|Modjaz et al.||2008| also used Mb values this facilitated 
easer adoption of and comparison with her data. 

3.3. Data Tables 

In Tables [T] [2] & [3] we present our line fluxes and esti- 
mated physicalgalaxy properties in tabular form. Line 
fluxes are presented with Galactic extinction removed 
but without correction for internal extinction. When the 
literature fluxes were not alr eady corrected f or Ga lactic 



extinction the value from the [Schlegel et al.| (|1998| dust 



maps was applied. We do this, since some metallicity 
determination methodologies also fit the inter nal extinc- 
tion simult aneously with the metallicity (i.e. Tremonti 
et al. 2004 1 , values are given without removing internal 



extinction. 

Internal extinction is calculated by determining an ex- 
tinction value which matches the measured Ha to H/3 
line flux ratio to the balmier decrement. If the Ha flux 
is not available but the H7 flux is, the extinction is sim- 
ilarly calculated from the H/3 to H7 ratio. The Ha flux 
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can then be estimated from the H/3 flux via the balmier 
decrement. This estimated Ha flux can then be used to 
obtain the SFR. 

In some cases, we have calibrated optical spectra that 
do not contain the Ha line, along with an uncalibrated 
near-infrared spectra of the Ha and 6583 A [N II] lines for 
the same object. In these cases, we can use the estimated 
Ha flux and the observed [N II] / Ha count ratio to 
estimate the 6583 A [N II] line flux. 

When the 4959 A [0 III] line was not available and the 
5007 A [O III] is, the 4959 A line is assumed to have | the 
5007 A lines flux as quantum mechanically required (see 
Graham et al.|2009[ ). Values given are measured without 
slit & fiber loss correction. The "spectrum" column de- 
notes whether we are presenting directly measured fluxes, 
fluxes normalized to H/3 or equivalent width values. If the 
H/3 column of an H/3 normalized flux is not 1 or 100 this 
is due to the galactic extinction correction being applied 
to the normalized values in the literature. The numeric 
superscript on the spectrum column indicated the refer- 
ence^) for the spectroscopy. Similarly, the superscripts 
on Mb, and redshift columns provide referen ces for those 
values. The reference key is given in Sect ion [3 . 3 . 1 1 which 
is separate from the general references of this paper. 

3.3.1. Tabular References 

1 Levesque, E. M., Berger, E., Kewley, L. J., & Bagley, 
M. M. 2010, AJ, 139, 694 

2 Svensson, K. M., Levan, A. J., Tanvir, N. R., Fruchter, 
A. S., & Strolger, L.-G. 2010, MNRAS, 405, 57 

3 Levesque, E. M., Kewley, L. J.. Berger, E., & Zahid, 
H. J. 2010, AJ, 140, 1557 

4 Sollerman, J., Ostlin, G., Fynbo, J. P. U., et al. 2005, 
New A, 11, 103 

5 Christensen, L., Vreeswijk, P. M., Sollerman, J., et al. 
2008, A&A, 490, 45 

6 Price, P. A., Kulkarni, S. R., Berger, E., et al. 2002, 
ApJ, 571, L121 

7 Garnavich, P. M., Stanek, K. Z., Wyrzykowski, L., et 
al. 2003, ApJ, 582, 924 

8 Hammer, F., Flores, H., Schaerer, D., et al. 2006, 
A&A, 454, 103 

9 Levesque, E. M., Kewley, L. J., Graham, J. F., & 
Fruchter, A. S. 2010, ApJ, 712, L26 

10 Levesque, E. M., Kewley, L. J., Berger, E., & Zahid, 
H. J. 2010, AJ, 140, 1557 

11 McGlynn, S., Clark, D. J., Dean, A. J., et al. 2007, 
A&A, 466, 895 

12 Graham et al. 2012 in prep 

13 Thone, C. C., Fynbo, J. P. U., Ostlin, G., et al. 2008, 
ApJ, 676, 1151 

14 computed from data given therein 

15 calculated from SDSS mags and redshift 

16 Modjaz, M., Kewley, L., Kirshner, R. P., et al. 2008, 
AJ, 135, 1136 

17 Ferguson, A. M. N., Gallagher, J. S., & Wyse, 
R. F. G. 1998, AJ, 116, 673 

18 Modjaz, M., Kewley, L., Bloom, J. S., et al. 2011, 
ApJ, 731, L4 

19 SDSS-mpg 

20 Sanders, N. E., Soderbe rg, A. M., Levesque, E. M., 
et al. 2012, |arXiv:1 206.2643 

21 NED 



22 Levesque, E. M., Soderberg, A. M., Foley, R. J., et 
al. 2010, ApJ, 709, L26 



4. ANALYSIS & RESULTS 

We begin our analysis by comparing the metallicity 
luminosity relation of LG RB and broad-lined Type Ic 
(Ic-bl) SNe hosts (S ection | 4.1[) as was originally done in 
the seminal work of Modjaz et al. (2008). Using a much 
smaller sample of LGRBs, |Modjaz et al. ( 2008[) provided 
some of the first strong evide nce that LGRBs show a 
strong metallicity aversion. 



|Modjaz et all (|2008l ) com- 
pared the luminosities metallicities ot LGRB ana Type 
Ic-bl hosts to the SDSS population. However, the SDSS 
suffers two major problems when used in this way as 
a comparison sample. First, it is not a volume limited 
sample, but a magnitude limited sample, and thus it is 
heavily weighted to bright (and therefore typically metal 
rich) galaxies. Secondly, it is primarily a very local sam- 
ple with median and maximum redshifts several times 
smaller than our LGRB sample, which extends nearly 
out to z ~ 1. We attack these problems with two ba- 
sic steps. Out of the main SDSS we create a volume 
limited survey, to allow comparis on wi th the (often very 
faint) LGRB hosts (see Section 
TKRS to provide a hig her redsr 
sample (see Section 2.2 1. Further 



2.1.11, and we use the 
shitt galaxy comparison 
we expand upon the 
Modjaz comparison sample of broad-line Type Ic SNe, 
by adding in Type II SNe, as well as a synthetic star- 
formation weighted population as further indicators of 
the distribution of star-formation. 

We do not rely upon the two-dimensional luminosity 
metallicity plot as the basis of our main analysis. In- 
stead, using the volume limited sample, we analyze how 
how each class of objects is distributed in metallicity with 
respect to its underlying star-formation and whether, 
once metallicity is accounted for, the probability of an 
object going off in a galaxy is proportional to the SFR 
of that galaxy. 

Our analysis shows that the offset LGRBs from the 
standard mass-metallicity relation for galaxies is robust. 
Inclusion of the bias towards lower metallicity caused 
by star-formation is not sufficient to explain the stark 
deviation of LGRBs. While the SNe track the star- 
formation distribution of the local universe, eighty per- 
cent of LGRBs are found in the ten percent of star- 
formation with the lowest metallicity. Nonetheless, the 
probability of an LGRB going off in a galaxy is directly 
proportional to its SFR, once metallicity has been taken 
into account. Indeed, due to an apparent similarity 
in star-formation distributions between low and higher 
metallicity galaxies, the distribution of LGRBs, just like 
the SNe, tracks the SFR of the entire sample of star- 
forming galaxies. Finally we show that the expected 
metallicity evolution with redshift of the LGRB host 
sample is much smaller then the LGRB metallicity aver- 
sion. Thus we conclude that the only remaining viable 
explanation for the metallicity distribution of LGRBs 
is an intrinsic preference for low (i.e. one-third solar) 
metallicity environments, which does not completely ex- 
clude LGRB formation at solar metallicity and above. 

4.1. Metallicity Luminosity Relation of Progenitor Sites 
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Table 1 

Table of LGRB hosts. Metallicity values are computed from the tabular line fluxes via the |Kobulnicky fc Kew lcy (2004) R23 diagnostic. 
When the 4959 A [O III] line is not listed and the 5007 A [O III] is available the 4959 A line is assumed to have | the 5007 A lines flux as 

quantum mechanically required (see |Graham et al.|2009"|f or a more detailed description). References for data in this table are indicated 
by numeric superscript with the key given m Section |a. 3. 1| and are separate from the general references of this paper. All given flux values 
have been corrected for galactic but not internal extinction (when the literature fluxes were not already corrected for galactic extinction 
the value from theJSchlcgcl ct al. 1998 dust maps was applied). The Spectrum column indicates type of spectral normalization with 

reference in superscript. Flux spectra are in units of 10 17 ergs s — 1 cm -2 
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Table 2 

Table of broad-lined Type Ic SNe hosts. The spectral values displayed are from sp ectra either of the the center of the galaxy or the site of 
the SN explosion. The metallicities are derived from these spectral values via the |Kobulnicky fc Kewley| < |2004] > R23 diagnostic when the 

required lines are available. Otherwise, the [N II] /Ha diagnostic is used. Metallicities directly calculated trom the tabular values are 
shown in bold. The metallicity values not shown in bold are shifted fro m the calculated (bold) values via the metallicity gradient method 
described in Appendix^] As in Table ^ references are given in Section 13 . 3 . 1 1 and are separate from the general references of this paper. 
Spectra column indicates type of spectra with reference in superscript. Flux spectra are in units of 10 17 ergs s — 1 cm -2 . The superscript 

by the SN name denotes whether the SN was detected in a targeted or non targeted manner as listed below: 
N SNe found in a non targeted manner 
T SNe found in a targeted manner 
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Table 3 

Table of Type II SNe hosts. Values are displayed as in Table [l] All of our Type II SNe have hosts in the SDSS Catalog. All of the values 

for these hosts, with the exception of the metallicity which is recomputed, come exclusively from the SDSS-mpg data products. Flux 
spectra are in units of 10 17 ergs s _1 cm -2 . The superscript a, b, & c by the SN name denotes the survey population in which the SN was 
originally detected as listed below (multiable listings indicate that the SN was detected in more then one population): As in Figure [2] the 

superscript "T" or "N" denotes whether the SN was detected in a targeted or non targeted manner respectively. 
a |Prieto et al.||2008| compilation (contains SNe found in both targeted and non targeted surveys) 
b SJJbS Supernova Survey (non targeted survey) 
c Nearby Supernova Factory (non targeted survey) 



As already noted, we begin our analysis by comparing 
the site metallicites and host luminosities of LGRBs and 
Type Ic-bl SNe with the central metallicities of the SDSS 
star-forming galaxy population. However, the SDSS is a 
very local sample. After the cuts based on galaxy emis- 
sion line signal to noise, described in Section |2.1.1| the 
median redshift of of our SDSS star-forming galaxy sam- 
ple is z = 0.0728, with no objects having a redshift above 
0.4 (with only few objects with a redshift above 0.3). Our 
LGRB population however, has a median redshift of z = 
0.36 and a maximum redshift of z = 0.83. Thus to extend 
the redshift range of the star-forming galaxy population, 
in this work we will also include the TKRS population 
described in Section 2.2 for much of our analysis. 



R23 metallicites are used for all LGRB hosts and star- 
forming galaxies, and SN hosts where available. For the 
fraction of SN hosts where R23 lines are not available, the 
less accurate [N II] /Ha metallicity is used. All metallic- 
ities were calculated by us using the ava ilable emission 
lines, and in a ll cases are placed on the Kobulnicky fc 
| Kewley| ( 2004 1 metallicity scale as described in Section 
|3.1| Due to the absence of line flux error data for many 
of the objects, a standard error of ±0.07 dex is assumed 
for R23 metallicities and ±0.15 f or [N II] /Ha metall icities 
was assumed. Compared to the Modjaz et al. ( 2008 ) sam- 
ple, our sample has an approximate three fold increase in 
the number of LGRBs, one additional SNe, and a second 
general star forming galaxy population (the TKRS) to 
better reflect the redshift distribution of the LGRBs. In 
order to avoid unknowingly biasing our SNe host prop- 
erties, we differentiate between SNe that were discovered 
in an untargeted search of the sky as opposed to SNe 



discovered in targeted surveys of specific galaxies. The 
properties of the hosts of the targeted SNe may be bi- 
ased by the selection criteria used in selecting the galaxies 
within which to search for SNe events. 

Nearly all LGRBs that have been studied sufficiently 
well to assure detection of an underlying Type Ic SN 
have shown evidence of a Type Ic-bl SN. Thus we might 
expect the progenitors of the general population of Type 
Ic-bl SNe to be closely related to those of LGRBs. A 
comparison of the hosts of LGRBs and those of Type Ic- 
bl SNe could give us some indication if this assumption 
is indeed true. 

In Figure [3] we plot metallicity vs. B band absolute 
galaxy luminosity for LGRBs, Type Ic-bl SNe hosts with- 
out associated LGRB events, and general star-forming 
galaxies (at low redshift fro m the SDSS and at h igh red- 
shift from the TKRS). As |Modjaz et al.| (|2008|) noted, 



the LGRB and SNe metallicity site values are more in- 
dicative of the progenitor local environment's metallicity 
than the galaxy core values and therefore we plot site 
values here, even though we only have core values for 
the background populations. We will address the issue 
of comparing between site metallicities of our LGRB and 
SNe populations and the central metallicites of our gen- 
eral star-forming galaxies in Section |4~2| 

Immediately obvious is the profound difference be- 
tween the LGRB metallicities and those of both the (non 
LGRB) broad-line Type Ic SNe a nd general star form - 

LET 



2008| 



sam- 



ing galaxies as originally shown by Modjaz e t al 
This effect is clear even though our larger LGRi 
pie contains three LGRB hosts at typical galaxy and 
SNe host metallicities. These objects, LGRBs 051022, 
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020819B, & 050826, are the focus of a companion pa- 
fl2012|) in prep. Nonetheless, the LGRB 



Graham 



per, 

population prefers much lower metallicity host environ- 
ments than the environments typically found in the gen- 
eral star-forming galaxy population or indeed the general 
Type Ic-bl SNe population. Thus, the primary conclu- 
sion that there is an intrinsic preference for low metal- 
licity in LGRB formation remains true. We will address 
this preference for low-metallicity statistically in Section 
IP1 



4.2. The Current Host Metallicity Luminosity Relation 

In this section we make two primary changes to the 
analysis already presented. We restrict the redshift range 
of the SDSS, and we add two new comparison popula- 
tions, Type II SNe and a star-formation weighted galaxy 
sample. Limiting the SDSS to the redshift cut of 0.02 < z 
< 0.04 described in Section [2. 1 . 1 1 avoids two major prob- 
lems. First, the SDSS is not a volume limited sample, 
but essentially a magnitude limited sample. It is there- 
fore heavily weighted to bright (and therefore typically 
metal rich) galaxies and thus neglects fainter galaxies 
more akin to LGRB hosts. Secondly the SDSS sample 
extends out to roughly z ~ 0.3, a range large enough 
that cosmic abundance evolution will start to become 
noticeable. However, with the narrow 0.02 < z < 0.04 
redshift constraint, the galaxies evolution issue is ren- 
dered negligible and the high completeness of the SDSS 
within this redshift range for galaxies brighter then -18 
converts the population to a volume limited survey 
allowing suitable comparison with the (often very faint) 
LGRB hosts. 

Unfortunately, sample size limitations preclude the 
universal application of this redshift cut across all popu- 
lations. Obviously the LGRB population cannot be use- 
fully constrained to the stated redshift range (this would 
reduce the LGRB hosts to a sample of one). Since evolu- 
tion of galaxies increases their metallicity over time this 
will introduce a low metallicity bias in the higher red- 
shift LGRB host population which we will subsequently 
demonstrate (using the TKRS sample) to be a much 
smaller effect than the low metallicity preference of the 
LGRBs in Section |4.4[ (Since the TKRS population is 



all at z > 0.3 they are not considered here). The broad- 
line Type Ic SNe sample also largely lies outside of the 
0.02 < z < 0.04 redshift range (with a significant fraction 
below). However since this population is naturally con- 
strained in redshift (median redshift of z = 0.0216 and a 
maximum redshift of z = 0.137), any biases in metallic- 
ity and host luminosity are limited, particularly as these 
objects have received deep follow-up observations. 

While the LGRB and Type Ic-bl SNe populations may 
themselves be (biased) tracers of star-formation (presum- 
ably a galaxy with twice the star- format ion would, all 
else remaining equal, have twice LGRB and SNe events) 
this is not the case for the SDSS population where the 
brighter galaxies have more star-formation but are still 
plotted singularly in Figure |J Thus, in Figure [IJ we also 
plot a sub-selection of the redshift cut SDSS population 
selecte d rand omly after weighting by star-formation (see 
Section 2.1.2 for details). 

Nature However, also provides a star-formation 
weighted galaxy sample through the hosts of Type II 
SNe. These SNe are probably largely unbiased with re- 



spect to metallicity (though their discovery will be af- 
fected by dust), and they are sufficiently abundant that 
we can obtain a reasonably-sized population of hosts 
which are entirely contained in the redshift-cut SDSS 
sample. This has the added benefit that the spectro- 
scopic observations and analysis are identical for the ob- 
ject hosts and the S PSS pop ulation. 

While in Section [0 and |Modjaz et al.| ( |2008[ ) the 
metallicity of the SNe and LGRB site within the host 
was used, here and subsequently we switch to using the 
central metallicities of the host galaxies to allow for a 
better comparison with the general star forming galaxy 
population. Since galaxies are not uniform metallicity 
objects, there will frequently be an offset in metallicity 
between the center of the host and that of the explosion 
site. In a few cases we have the metallicity of both the 
explosion site and the center of the host. In most cases 
do not have a measurement of the metallicity offset, but 
can estimate it using the method described Appendix [A"] 
. Examining Table [2J one sees observed and expected orN 
sets up to 0.3 dec in metallicity between the explosion 
site and the host center fo r the broad-line Type Ic SNe. 

As discussed in Sec tion |2.4[ LGRBs ten d to occur in 
the brightest regions (Fruchter et al. 2006) of their pre- 
dominately irregular host galaxies. Furthermore due to 
the distance to the LGRB hosts, and their generally com- 
pact sizes, most of the light of the ho st goes down a single 
ground-based slit (see Section 3.2.1| . Thus in most cases 



we would not expect a significant offset in metallicity 
between the host site and the measured host metallicity. 
However, in the case of the three well-resolved LGRB 
hosts, LGRBs 980425, 020819B, & 060505, we are able 
to directly address this issue, and we discuss their (small) 
metallicity offsets in Section [2~4| 

In examining Figure [4j one notes that the agreement 
between the luminosity - metallicity relations of the SNe 
host and the red-shift cut star-forming SDSS sample is 
notably better than the agreement with the entire SDSS 
star- forming population (shown in Figure [3]). The agree- 
ment of the two-dimensional distribution is even bet- 
ter when one compares the SNe to the star-formation 
weighted sample (small blue dots in Figure [4|. The SNe 
and the synthetic star-forming galaxy population appear 
to follow the same luminosity-metallicity relation and 
have similar distributions of population density in the 
two-dimensional luminosity metallicity plane. Nonethe- 
less, the LGRBs remain well-separated from all of the 
comparison galaxy samples. In the following section, we 
will use a more powerful statistical methodology to com- 
pare the properties of the hosts of the SNe, LGRBs, and 
the star-forming SDSS galaxy population. 

4.3. A Statistical Comparison of Populations 

In the previous section, we examined the host metallic- 
ity distribution with respect to luminosity through the 
use of the metallicity luminosity plot. Here we take a 
more statistical view and ask specifically: do LGRBs and 
SNe track star-formation independent of host metallic- 
ity? In the second section we have a first look at the 
question: does the host SFR itself bias the probability 
of LGRB or SNe star-formation, as one might expect if 
star-burst conditions enhanced the production LGRBs 
or SNe beyond the SFR itself, perhaps by altering the 
IMF towards higher mass objects? 
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Figure 4. Mctallicity vs. B band absolute galaxy luminosity. The Type II SNe hosts (purple triangles), SDSS galaxies (small green 
points) and the star-formation weighted random SDSS galaxi es syn thetic population (small blue points) shown here are limited to the 
sample redshift range of 0.02 < z < 0.04 discussed in Section |2.1.1| In order to get reasonable comparison samples, no redshift cut has 
been applied to the Type Ic-bl host population (blue circles) or the LGRB population (red squares). This may produce a small bias toward 
brighter hosts for the broad-line Type Ic SNe; however as the most distant broad-line Type Ic SNe has a redshift of 0.137, it will not produce 
a significant bias in metallicity. Fo r LG RBs, to some extent, both biases will occur and we discuss the small effect of an unrestricted redshift 
distribution on LGRBs in Section [O] As in Figure [3] filled and open symbols represent SNe selected in untargeted and targeted surveys 
respectively. In ord er to e xplicitly trace star-formation we introduce the star-formation weighted synthetic population (small blue points) 
discussed in Section |2.1.2| This population tracks the SNe quite well in luminosity-metallicity space, as one would expect with SNe tracking 
star-formation. 
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4.3.1. Star-formation Metalhcity Distribution 

If LGRBs and SNe events track star-formation, then 
their distribution in metalhcity space should be identical 
to that of star-formation in general. Thus if in the local 
universe twice as many stars are generated per unit time 
at metalhcity Z\ than Zo, then we should hnd twice as 
many SNe and LGRBs at a metalhcity of Z\ than at Zq. 
To look at this issue we think of each host galaxy, or 
galaxy in the SDSS, as a particle of given metalhcity and 
star-formation. We then sort our objects, and galaxies by 
metalhcity, and examine what fraction of star-formation, 
or object formation, occurs up to a given metalhcity. 

For the SDSS, we create a normalized cumulative dis- 
tribution of star-formation versus metalhcity, i.e. sort 
galaxies by their metalhcity, and determine the fraction 
of total star- format ion in the SDSS population contained 
in galaxies below a given metalhcity. This cumulative 
distribution can be seen in Figure [5] For the LGRBs 
and SNe we assume the null hypothesis that these ob- 
jects track star-formation, and thus each explosion can be 
thought of as a "quantum" of star- format ion. In Figure 
[5] we can therefore plot their cumulative distributions as 
a step function of objects. A sim ilar LGRB step function 



( |2006[) to compare the 
cm of he 



was hrst used in |Fruchter et al 

burst location with the distribution of host galaxy light, 
and has more recentl y been applied to the me talhcity of 
SNe independently by |Kelly fe Kirshner| ( |2012| ). The pri- 
mary innovation here is the direct comparison between 
the step functions of events with a normalized cumulative 
distribution function of the star- format ion in the general 
galaxy population. This analysis has the notable advan- 
tage of allowing comparison of events assumed to track 
star-formation without actually requiring information on 
the SFR of the hosts of these events themselves. In our 
case this is primarily of benefit with the Type Ic-bl SNe 
and the LGRB's where such information is not readily 
available for a number of the objects. 

The general agreement of both the Type Ic-bl and Type 
II SNe with the SDSS general star-forming galaxy popu- 
lation seen in the plots is what one would expect assum- 
ing the SNe are a product of the star- format ion alone 
without additional selective pressures. Indeed, if we re- 
strict our SNe samples to only hosts with Mb brighter 
than -18 magnitude, which is the completeness limit of 
our SDSS sample, all of the SNe populations agree rea- 
sonably well with the SDSS distribution of star-formation 
with KS values greater than ^0.1 in all cases. It is widely 
realized that targeted surveys will have selection effects. 
Although our nearby sample may somewhat limit this 
effect, in general, targeted surveys overrepresent bright 
(and attractive) galaxies. However, even untargeted sur- 
veys can have significant biases, and those used to cre- 
ate our non-targeted SNe samples are no exception. For 
example, the Palomar Transient Factory (PTF), which 
contributes some of our Type Ic-bl SNe, has placed sig- 
nificant emphasis on typing SNe found in faint or missing 
hosts (Y. Arcavi private communication). This will bias 
metallicities due to the mass-metallicity relationship of 
galaxies. But perhaps the greatest issue for our untar- 
geted samples is due to the difficulty of finding SNe on 
top of bright surface brightness hosts, or near the nuclei 
of massive hosts. We find that this effect can be charac- 
terized fairly well by examining the fiber plug magnitudes 



from the centers of the SDSS galaxies, which act as an 
estimate of the galaxies' central surface brightnesses (see 
Appendix Ib]) . Interestingly, the SNe sample that most 
closely tracks the SDSS star-formation metalhcity func- 
tion (KS — 0.88) is the combined sample of untargeted 
and targeted Type Ic-bl SN in hosts brighter than -18 
(see Figure [5] right plot). 

However, the LGRB population differs greatly from 
the distribution expected based on star-formation. In- 
deed a K-S analysis suggests a probability of less than 
10~ 5 ! From this comparison it is apparent that LGRB's 
occur much more preferentially than star-formation in 
low metalhcity environments. As can be seen in Figure 
[5j while untargeted and targeted samples have differing 
biases (we will discuss possible biases in the untargeted 
population in greater detail in the next section), those 
biases are nowhere near sufficient to reproduce the devi- 
ation from star-formation seen in the LGRB population. 
Indeed, 11 out of 14 LGRBs lie in the lowest metalhcity 
tenth of the star-formation in the volume-limited SDSS 
sample. This is in spite of the fact that the LGRB sample 
is not volume limited, and thus the requirement that we 
be able to measure the metalhcity of the host, biases us 
towards including more luminous, and thus likely more 
metal rich, galaxies. There is one bias on the LGRB pop- 
ulation that tends toward lower-metallicites. A fraction 
of our LGRB sample, however, extends up to a redshift 
of one. This could bia s us to lower metallicites. However, 
as we show in Section [4~4} this bias is not at all sufficient 
to produce the observed effect. 

4.3.2. Fractional star-formation analysis 

Here we employ an alternate analysis directly compar- 
ing the star-formation of target host galaxies with those 
of the SDSS general star-forming galaxy population. We 
ask for each SN or LGRB, what fraction of star-formation 
in the SDSS is contained in galaxies of equal or lesser 
star- format ion than the host. We do this to determine 
whether there is any apparent bias towards highly star- 
forming galaxies, beyond a linear correlation with SFR. 

We compare each of the hosts both to the SDSS pop- 
ulation, and also to a subset of the SDSS galaxies with 
comparable metallicities. However, we cannot just use 
the volume-limited redshift cut SDSS population for this 
procedure. The LGRB sample is essentially a magni- 
tude limited sample and is thus biased towards bright, 
and therefore highly star-forming galaxies. The entire 
SDSS star-forming galaxy sample, however, is a magni- 
tude limited sample and thus is the correct comparison 
with the LGRBs. The Type II SNe are volume limited 
to the same range as the SDSS redshift cut and thus the 
volume limited SDSS sample is applied. We also include 
in our comparison the Type Ic-bl SNe from our sample 
which have SDSS spectroscopy. This allows us to obtain 
a SFR for their hosts. However, we make no redshift cut 
on this sample, and allow redshifts both below and above 
the limits of the redshift cut SDSS sample. Nonetheless, 
the Type Ic-bl SNe are far more magnitude limited by 
the brightness of their SNe than they are by their host, 
and thus the volume limited SDSS is a reasonable, but 
imperfect, comparison sample. We therefore adopt it for 
determining the Type Ic-bl fractional star-formation val- 
ues. 

In Figure [6] we show the results of this procedure. In 
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Figure 5. Cumulative fraction of population (LGRBs and SNe) or total star-formation (redshift-cut SDSS star-forming sample) vs. galaxy 
central metallicity. LGRB host galaxies are shown in red, Type Ic-bl SNe hosts in blue, and Type II SNe hosts in purple. Non-targeted 
SNe are shown with a dashed line, targeted SNe with a dotted line and the union targeted and non-targeted populations with a solid line. 
The redshift-cut SDSS star-forming galaxy population is shown in black. Only the SDSS galaxies, and the Type II SNe hosts (which are all 
included in the SDSSS sample), are limited to the redshift cut range. In the left-hand plot we use our full samples (subject to the redshift 
cuts). In the right-hand plot we restrict both the host and SDSS samples to the magnitude the approximate completeness limit of the 
SDSS in our redshift range of Mg < —18 to insure that host magnitude does not bias the metallicity selection. A detailed description of 
the plotted populations is described in Section|2] The LGRB hosts have metallicities considerably lower then would be obtained simply by 
following star-formation while both types of supernova are consistent with the star-formation distribution. Thus the metallicity distribution 
of LGRBs cannot be explained only by association with star-forming galaxies. 



the left-hand plot, for each object (LGRB, SNe, etc.) we 
divide the total star-formation in SDSS galaxies with in- 
dividual SFRs less then or equal to that of the target's 
host galaxy with the total star-formation in SDSS galax- 
ies to create a fractional star-formation value. In the 
right-hand plot, the same procedure is used, but hosts are 
compared only to galaxies of similar metallicity, rather 
than the entire star-forming sample. These sorted frac- 
tional values are shown in a standard normalized cumu- 
lative distribution plot. A diagonal reference line corre- 
sponding to perfect tracking of the star- format ion dis- 
tribution is shown for comparison. One would expect a 
distribution of objects that followed star-formation irre- 
spective of metallicity to agree well with the diagonal of 
the left-hand plot. However, a deviation from that line, 
might be expected of a population that followed star- 
formation, but was also biased by metallicity. However, 
such a population should follow the diagonal line in the 
left-hand plot, as should a population, that has no metal- 
licity bias. The synthetic star-formation population was 
tested against the diagonal reference line as a check on 
methodology and implementation and was found to be 
is good agreement, but we do not show the synthetic 
population in this plot. 

With the exception of the untargeted Type II SNe, 
all of the SN and LGRB population shown reasonable 
agreement with the general star-forming galaxy popula- 
tion (when the latter population is weighted by SFR) 



shows that the individual SFRs of these host galaxies 
tracks the general population. Their distributions are 
representative of events occurring as would be expected 
where they simply a product of star-formation. For the 
non-LGRB SNe this appears to be the entirety of the 
explanation necessary, however the trend of LGRBs to 
occur in low metallicity hosts means that while metal- 
licity is obviously a critical factor in burst occurrence 
the star-formation distribution does not alter nor is al- 
tered by the low metallicity nature of these host galax- 
ies. Thus the observed metallicity bias can not be the 
incidental product of a deviation towards a host galaxy 
population of atypical SFR. The untargeted Type II SNe 
appear to suffer from a strong bias against being found 
on high-surface brightness galaxies, and we believe this 
is the cause of their deviation from the expected star- 
formation track. The evidence for this bias is discussed 
in detail in Appendix |B} 

As a further check to insure that metallicity has not 
biased our results, in the right hand plot of Figure [6j 
we isolate the underlying star-formation by removing 
the metallicity component from the analysis. To effect 
this we reduce the SDSS sample used for the fraction 
value computation to only those galaxies whose metal- 
licity matches (± 0.1 dex) of the target object. Thus for 
a low metallicity LGRB we would limit comparison to 
other low metallicity objects. This isometallicity analy- 
sis allows a comparison of the star- format ion tracking as 
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Figure 6. Normalized cumulative distribution of star-formation fractional values for the SNe and LGRBs. (see Section |4. 3. 2\ . For each 
object we determine the fraction of star-formation contained in SDSS galaxies of equal or lesser star-formation than the object's host. 
Were the objects distributed as star-formation the objects would be expected to follow a diagonal line on these plots. The left plot uses 
star-formation across all mctallicities to determine the fractional values whereas in the right plot uses only star-formation within ±0.1 dex 
of the metallicity of the object. For the SNe we use the volume limited SDSS sample as a comparison. However, as the LGRB host sample 
is itself approximately a magnitude limited sample, the LGRBs are compared to the total SDSS star-forming sample. These comparison 
choices are discussed in more detail in the main text. The targeted Type II SNe (blue line) and LGRBs (red line) track the diagonal quite 
well indicating that both populations have fairly typical star-formation distributions. However a similar comparison of the non-targeted 
Type II SNe sample (solid green line) and the Type Ic-bl SNe (grey line) show lower then expected host SFRs. This is especially apparent 
when considering only Type II SNe found in the SDSS and SNFactory searches (dashed green line). We suspect that this is due to 
limitations in the ability of non-targeted searches to identify SNe events on high surface brightn ess backgrounds, as is d iscussed in greater 
detail in Appendix Jb] Thus while LGRB hosts do have higher than average SFRs, as noted by |Mannucci et al.|l|2011| , their SFRs agree 
with what one woufa expect of a population that tracks star-formation. Thus metallicity, rather tnan BFrl, must be tne 
the discordant results shown in Figure [5] 



: the primary source of 



influenced by all remaining factors. The good agreement 
between the LGRB, targeted Type II SNe population, 
and the general star-forming galaxy population (when 
the latter population is weighted by SFR) shows that 
the individual SFRs of these host galaxies are typical of 
galaxies of similar metallicity, and indeed of the entire 
population, as shown previously. 

The original similarity between the star-formation dis- 
tribution of LGRB hosts, a population heavily biased 
towards low metallicity, and SDSS galaxies, whose star- 
formation we showed in the previous section to be pri- 
marily at high metallicity is perhaps surprising. We fur- 
ther explore this by plotting the SFR distribution of the 
low metallicity SDSS galaxies weighted by their star- 
formation against the similarly weighted star-formation 
distribution of the entire SDSS sample (see Figure [7]). 
This is equivalent to plotting fractional values (just as 
in Figure p] left) for a low-metallicity subset of the star- 
formation weighted galaxy sample. From this we find 
that the SFR distribution of the low-metallicity SDSS 
galaxies is reasonably similar to that of the entire SDSS 
population. This is actually done twice to show that the 
similarity is present in both our full SDSS sample and the 
volume limited subset. We also overplot the LGRB pop- 
ulation on FigureJTJshowing the full SDSS sample. This is 
done with both the full LGRB sample and after discard- 



ing the three high metallicity objects to better match the 
low metallicity SDSS population. However the LGRB 
population is too small to differentiate between the star- 
formation distribution of the low-metallicity SDSS galax- 
ies or the full SDSS galaxy sample. 

4.4. Host Metallicity Luminosity Redshift Relation 

While our SNe samples are all restricted to local events, 
LGRBs tend to be cosmological, and thus local events are 
rare, precluding a similar redshift restriction. Thus our 
LGRB sample extends out to nearly a redshift of one, 
while our volume limited SDSS sample only extends to 
z < 0.04 and the entire SDSS to about z < 0.4. While the 
intrinsic metallicity preferences of LGRB and SNe should 
be independent of redshift, the metallicities of galaxies in 
general change with redshift, and thus any understand- 
ing of the effects of metallicity on LGRB formation must 
also include the possible biases introduced by the metal- 
licity evolution of galaxies with redshift. To extend our 
general star-forming galaxy comparison population out 
to a comparable redshift as the LGRB population we 
therefore add the TKR S galaxy survey population as de- 
scribed in Section [272] The TKRS population covers the 
redshift range 0.3 < z < 1.0, and thus the the SDSS and 
the TKRS together span the entire redshift range of the 
LGRB sample. 

In Figure [8j we plot the host central metallicity lu- 
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Sorted Fraction of Total Population Star-Formation Sorted Fraction of Total Population Star-Formation 

Figure 7. Normalized cumulative distribution of star-formation of low metallicity galaxies (log(0/H)+12 < 8.4) versus the entire SDSS 
population. Left plot is for the volume limited SDSS population, right plot is for the entire SDSS population with the LGRB population 
(red lines) overplotted. For each low metallicity SDSS galaxy we determine the fraction of star-formation contained in SDSS galaxies of 
equal or lesser star-formation in the low metallicity sample and in the entire SDSS. We then plot the former as the abscissa and the latter 
as the ordinate. The two samples are perhaps surprisingly similar, though as one might expect, the larger SDSS sample has a somewhat 
higher fraction of extreme star-forming galaxies. Both the entire LGRB host sample and the subset of low- metallicity hosts are both plotted 
(solid vs. dashed red line respectively). The star- formation distributions are sufficiently similar (and our LGRB sample sufficiently small) 
that neither is obviously a better match for the LGRB population. 



minosity redshift relation in a similar manner to previ- 
ous plots, but instead of using color to separate object 
classes, we instead us it to index the redshifts of the ob- 
jects. In our previous plots, we implicitly ignored the 
effect of galaxy chemical evolution over time, and thus 
assumed this effect to be much smaller than the metallic- 
ity offset of LGRBs from the general star-forming galaxy 
population. Here we explicitly show this to be the case. 

The different redshift distributions of the SNe and 
LGRBs populations mean that they are best compared 
not with each other but with the general star-forming 
pop ulations at similar redshifts. As shown in Section 
4.2 the SNe are consistent with the general star-forming 
galaxy population at similar redshifts. The LGRB pop- 
ulation however is not consistent with the general star- 
forming galaxy population at any redshift. Despite the 
existence of three relatively high metallicity LGRB hosts, 
the LGRB population is on average at far lower metallic- 
ities than the distribution of our comparison galaxy sam- 
ples across any redshift range It should again be noted 
that the TRKS population above a redshift of z ^ 0.5 
(where the [N II] and Ha lines move into the infrared) 



is assumed by |Kobulnicky fc Kewley 



( |2004[ ) to be on the 
upper branch of the R23 metallicity diagnostic. Thus the 
~ 10% of the sample that would b e exp ected to be on the 
lower branch is missed (see Section 2.2 for a more detailed 
discussion of this issue). However it is apparent in Fig- 
ure [8] that a ten percent correction in the TKRS galaxy 
numbers is insufficient the bring the galaxy population 
anywhere near agreement with the metallicity distribu- 
tion of the LGRB hosts. 



Interestingly, the three high metallicity LGRBs 
(051022, 020819B, & 050826) are consistent with the 
general star-forming galaxy population for galaxies of 
their luminosity and redshift. If the metal aversion effect 
was evident in the high metallicity LGRB population we 
might expect them to be close to, but still biased below, 
comparable galaxies. However, for the three cases avail- 
able they seem to be neither outliers with respect to the 
comparison population or even among the lowest galaxies 
available within them. It therefore seems possible that 
above some metallicity, any further bias towards lower 
metallicities is muted. 

5. SUMMARY 

In this paper, we compare the metallicity distribu- 
tion of the hosts of Long duration Gamma Ray Bursts 
(LGRB) with that of several other related populations: 
Type II supernovae (SNe) (which can serve as markers of 
star- formation) , broad-lined Type Ic (Ic-bl) SNe (which 
as discussed earlier are closely related to LGRBs), and 
the general star-forming galaxy population. We perform 
this comparison across these different populations not 
primarily by number of galaxies but rather by weight- 
ing by their star- format ion rates (SFRs). Our results 
show that, not only are the hosts of LGRBs at lower 
metallicites than either SNe hosts or general star-forming 
galaxies, but also that while Type Ic-bl and Type II 
SNe track star- format ion (within our statistical ability 
to measure), more then three-quarters of our LGRBs 
are clustered in the lowest metallicity tenth of the star- 
formation. 
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Figure 8. Central metallicity vs. B band absolute magnitude of LGRB (squares) and Type Ic-bl SNe hosts (circles). As before (figures 
\3\ and HI filled & open circles represent SNe selected in an untargeted and targeted manner respectively. Star-forming galaxies from the 
&DSS (small dots) and TKRS (diamonds) are plotted in the background to provide low and high redshift comparison samples respectively. 
While in the previous plots, objects are color coded by type, here color is used to index redshift. The different redshift distributions of 
the SNe and LGRBs populations mean that they are best compared not with each other but with the general star-forming populations at 
similar redshifts. While as discussed earlier, the SNe we'll match their comparison sample; however, with the exception of a few relatively 
high metallicity LGRB hosts, which match the metallicities of similar magnitude galaxies in the TRKS, the LGRB population as a whole 
remains at far lower metallicities than even the TKRS. 
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To the end, we have complied all spectra for LGRB and 
Type Ic-bl SNe hosts with host emission spectroscopy 
sufficient to allow metallicity measurement. To relate 
these events to the general star-forming galaxy popu- 
lation we have extracted the approximately 137 thou- 
sand galaxies from the Sloan Digital Sky Survey (SDSS) 
with line measurements for suitable metallicity measure- 
ment. This has also allowed us to compile a sample of the 
more common Type II SNe hosts simply by selecting the 
galaxies within our existing volume limited SDSS sample 
which are known to have hosted such Type II SNe events 
with all the expected advantages of inter-sample consis- 
tency. Additionally, in order to better match the redshift 
range of our LGRB population, we extend our general 
star-forming galaxy population out to a redshift range of 
approximately unity via the addition of the higher red- 
shift Team Keck Redshift Survey (TKRS) galaxy popu- 
lation. 

Our analyzes are based on four physical properties 
redshift, metallicity, rest frame B band absolute mag- 
nitude, and (for most of our galaxies) the SFR. To max- 
imize inter-sample consistency we have calculated inde- 
pendent metallicities using a common metallicity diag- 
nostic, scale, and code via the R23 method. The method 
uses the ratio of oxygen to hydrogen line strengths (and 
the ratio of two oxygen lines to characterize the degree 
of ionization) to estimate the oxygen abundance in HII 
regions as a proxy of the total galaxy metallicity. R23 
however has a metallicity degeneracy issue where metal 
line emission cools the election temperature causing the 
oxygen line strength, that was originally increasing with 
metallicity, to subsequently decrease yielding undiffer- 
entiable oxygen line strengths for both a high and low 
metallicity value. To resolve this degeneracy we have 
used observations of the [N II] /Ha line ratio, itself a 
crude metallicity indicator, to select between the upper 
and lower metallicity branch. For a few Type Ic-bl SNe 
we are forced to rely on this [N II] /Ha diagnostic ex- 
clusively with a corresponding increase in error. When 
comparing galaxies we use the B band galaxy luminosity 
rather than galaxy mass as reliable host mass estimates 
were not available for a substantial fraction of the ob- 
jects in our samples. We dete rmine the S FRs f rom the 
galaxy's Ha emission via the Kennicutt ( 1998 1 metric. 
The primary difficulty with this is adjusting the Ha line 
flux for slit and fiber losses. 



As shown in Modjaz et al. (2008) the low metallicity 
bias of LGRBs is visually apparent in a simple scatter- 
plot of host metallicity versus absolute magnitude, when 
compared to that of the hosts of Type Ic-bl SNe and 
the general star-forming galaxy population, represented 
by the SDSS. This result becomes even more impressive 
with the approximate three fold increase the the num- 
ber of LGRBs and a 50 % increase in number of Type 
Ic-bl SNe presented in Figure [3j In this figure, we also 
add a second general star forming galaxy population, the 
TKRS, to better reflect the higher redshift distribution 
of the LGRBs. This point is discu ssed in much greater 
detail in Figure [8] and Section |4~4") In our first presenta- 
tion of the scatter plot, we present the metallicity mea- 
sured at the location of the LGRBs and SNe within their 
hosts to show the actual metallicity of the environment 
responsible for their creation. In later plots, however, 
we switch to plotting the host galaxy central metallici- 



ties in order to better compare with the general galaxy 
population. To obtain a fair sample of the galaxy pop- 
ulation, we use a subset of the SDSS star-forming pop- 
ulati on res tricted to a volume limited sample (see Sec- 
tion 2.1.1 for details). Even this population, however, 
does not represent the way we expect SNe or LGRBs 
to choose their hosts. Their probability of going off in 
a particular galaxy should be proportional to the rate 
of star-formation in that galaxy (all other things being 
equal). To better emulate the expected occurrence of 
LGRBs and SNe, we then further select among these 
galaxies via random selection weighted by the underlay- 
ing star-formation of each galaxy. These volume limited 
and star-formation weighted samples are plotted along 
with the similarly volume limited Type II SNe in Figure 

m 

As can be seen in Figure HJ the distribution of the 
star- format ion weighted sample is significantly different 
from the distribution of galaxies by number with a large 
bias towards brighter galaxies having a greater fraction 
of star-formation. As the LGRBs in our sample extend 
to a redshift of about one, far beyond the reach of the 
SDSS, let alone our smaller volume limited survey, we 
must consider the evolution of metallicity with redshift 
as a complicating factor. Thus, later in the paper, in 
Figure [8l we expand our plotting by indexing object red- 
shift with color to provide a visual representation of the 
evolution of metallicity over time. From this plot it is 
clearly apparent that the effect of abundance evolution 
over the redshift range of our sample is nowhere near as 
dramatic as the metal aversion of the LGRB population 
as a whole. While there are three LGRB hosts in the 
high metallicity range of our sample (LGRBs 051022, 
020819B, & 050826), is it clearly apparent in all three 
metallicity vs. luminosity scatterplots (figures [3j [4j &[8| 
that while the SNe are distributed with the general star- 
forming galaxy population the LGRBs are on whole at 
the bottom of the metallicity distribution. It is notable 
that the three high metallicity (out of 14 total) LGRBs 
do appear consistent with the general star forming galaxy 
population of comparable brightness and redshift. This 
is intrinsically surprising as were the metal aversion to 
remain in effect for these objects we would expect their 
occurrence, if still in the high metallicity range, to be far 
lower then the typical metallicity for the population at 
that luminosity and redshift (i.e. either a outlier of said 
population or among the lowest galaxies available within 
it). 

While the metallicity vs. luminosity scatterplots show 
quite conclusively that the hosts of LGRBs are at lower 
metallicites than either SNe hosts or general star-forming 
galaxies, they do not directly address whether this metal- 
licity bias could be caused b y the anti-correlation be- 
tween SF R and metallicity of Mannucci et al. (12010 ) as 

iirc trT 



claimed in Mannucci et al. (20 



Mannucci et 
l]) . To contro 



is issue, 



we plot in Figure [5] tne integrated star- formation of the 
SDSS sample as a function of metallicity in comparison 
with the cumulative distributions of LGRBs and SNe. A 
similar plot, with a much smaller sample LGRB s ample 
and no com parison of SNE, was first shown by Stanek| 



etaL| (|2006| ). Both the non LGRB Type Ic-bl SNe and 
Type ITSNe track the distribution of star-forming SDSS 
galaxies quite well. The LGRBs however display a pro- 
found preference for lower metallicities. Thus the CCSNe 
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track star-formation independent of its metallicity, while 
LGRBs do not. 
However, this result could still be consistent with the 



Mannucci et al | (|2010j relation, if the SFRs of the LGRB 
host were wildly discrepant from the other populations. 



6\ (left plot) we di- 
and SNe hosts to 



To exclude this possibility, in Figure 
rectly compare the SFRs of the LGRi 
the general star-forming SDSS galaxy population. This 
comparison is performed by taking the SFR of each 
LGRB and SNe host and determining the fraction of the 
total star-formation in the general SDSS galaxy popula- 
tion that occurs in galaxies with less star- format ion than 
the host. These fractional values are then sorted and 
plotted as a normalized cumulative histogram. Should 
the distribution of star-formation for a given object type 
follow the general star-forming SDSS galaxy population 
then this histogram would track a diagonal line on the 
plot. For the SNe, the SDSS comparison sample is vol- 
ume limited. However, for the LGRBs, which are in- 
trinsically a magnitude limited sample, we use the entire 
SDSS, as a comparison magnitude limited sample. 

Both the LGRBs and targeted Type II SNe population 
track the diagonal well, indicating a good correspondence 
between the SFRs of the two populations and the gen- 
eral SDSS galaxy population, and suggesting that SFR 
is correlated with LGRB and SNe formation. However, 
as further refinement we can also generate fractional val- 
ues for the LGRB and SNe hosts using only a subset of 
the SDSS galaxy population with similar metallicity to 
each host whose fractional value is being calculated. This 
similar metallicity comparison, shown in Figure (right 
plot) , yields a slightly closer agreement between the SFR 
distribution of the LGRBs and the SDSS galaxy popula- 
tion. 

Given the low metallicity nature of the LGRB popula- 
tion the absence of any profound differences in these dis- 
tributions is itself interesting. We explore this in Figure 
[7] by plotting the SFR distribution of the low metallic- 
ity star-formation versus the star-formation of the entire 
population. One way to show this, would be to plot a 
low-metallicity subset of the randomly generated star- 
formation weighted galaxy population just as we have 
done for the LGRB and SNe hosts. However, this then 
relics on our generating a very large number random 
population. We can achieve the same effect, without 
relying on random selection, by using the whole sub- 
set of low-metallicity star-forming galaxies and compar- 
ing their fractional star-formation values in that subset 
against their values in the whole population, which is 
what we actually do in the figure. From this we find 
that the SFR distribution of the low-metallicity SDSS 
galaxies is reasonably similar to that of the entire SDSS 
population and because of this similarity , the LGRB 
population is too small to determine whether it prefer- 
entially tracks the SFR distribution of the low-metallicity 
SDSS galaxies or the full SDSS galaxy sample. Curiously 
the untargeted Type II SNe population fails to track the 
expected diagonal, even more so when the untargeted 
population is restricted to objects found trough untar- 
geted SNe surveys (as apposed to incidental detections 
from other untargeted observations). We believe this bias 
to be an issue with at least one of the untargeted sur- 
veys failing to detect SNe in high surface brightness back- 
grounds (see Appendix [5] and Figure 11). 



Thus the populat ion distribution suppor ts the scatter 
plots and rules out Mannucci et al. (20101 as an expla- 
nation. 

6. CONCLUSIONS 

In this paper, we have shown quite conclusively in var- 
ious scatterplots that the hosts of LGRBs are at lower 
metallicites than either SNe hosts or general star-forming 
galaxies. How ever, in order to directly address the ques- 
tion posed by Mannucci et al. (20111), of whether SFR 



can explain metallicity, we plot as function of metallicity 
the integrated star-formation of the SDSS sample in com- 
parison with the cumulative distributions of LGRBs and 
SNe. There we find that three quarters of our LGRBs 
are found in the in the bottom low-metallicity tenth of 
the star-formation, with the remaining quarter (like our 
SNe populations) appearing to track star- format ion in- 
dependent of metallicity. Assuming, all else being equal, 
that LGRBs (and SNe) occur proportionally to the al- 
lowable star- format ion this bias indicates that LGRBs 
clearly prefer much lower metallicity ho st environments 
(as first suggested in |Stanek et al. ( 2006 ) and shown with 
a much larger sample here). We also consider the SFR 
distribution of the LGRB and SNe populations and find 
that they are consistent with the SFR distribution of the 
SDSS sample and especially consistent with this compar- 
ison is limited only to SDSS objects of similar metallic- 
ity. Interestingly we find that the SFR distribution of 
the low-metallicity SDSS galaxies is reasonably similar 
to that of the entire SDSS population and cannot deter- 
mine which the LGRB population best tracks. 

In contrast with LGRBs, both our Type II and Type 
Ic-bl SNe populations appear to track the metallicity 
of the integrated SDSS star-formation. This is what 
one would expect if star- format ion alone is sufficient to 
explain the metallicity distribution of both SNe types. 
Nonetheless, any metallicity bias, even one much more 
minor than that seen with LGRBs, has significant impli- 
cations for a unified standard formati on model for Type 
Ic SNe Type Ic-bl SNe and LGRBs. [Kelly fc Kirshner 
( |2012[ ); jSanders et al.l 12012| ) have looked at the metal- 
licity distribution ot airterent SNe types vs. each other 
(without relation to the SDSS or LGRBs) and highlight a 
bias towards the Type Ic-bl SNe preferentially occurring 
in slightly lower metallicities than the other SNe popu- 
lations (including the non Type Ic-bl SNe). Due to our 
strict volume limits, imposed to allow comparison with 
the SDSS, our Type II SNe sample is far smaller than 
that used in these other works. And we have made a 
decision to only look at two SNe types — Type II SNe 
as hopefully a truly unbiased indicator of star- format ion, 
and Type Ic-bl SNe, due to their close association with 
LGRBs. Thus, while our work is well suited for its 
primary purpose of determining whether the metallic- 
ity distribution of LGRB hosts differs from other star- 
forming galaxies, it is much less powerful for distinguish- 
ing relative metallicities preference of SNe. However, the 
metallicity differences between Type Ic-bl and other SNe, 
hinted at in our sample and perhaps seen more clearly in 
these other works, is dwarfed by the strong metallicity 
bias seen in LGRBs. 

The presence of a metallicity bias between LGRBs and 
Type Ic-bl SNe poses a problem for explanations of the 
LGRB metallicity bias being the incidental result of an 
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IMF difference in their host galaxies. In nearly all cases 
where one would have expected to detect a Type Ic un- 
derlying a LGRB to have been detected, one has been 
found, and where good sp ectros c opy i s available, the 
Type Ic is broad-lined (c.f. Cano (2012) for a good dis- 
cussion of this point). Thus, many and perhaps all Type 
Ic-bl SNe share a common progenitor type with LGRBs. 
One would therefore expect the masses of Type Ic-bl SNe 
progenitors to be similar to those of LGRBs. Thus while 
IMF differences between galaxies could play some role in 
determining where one finds LGRBs, the IMF is almost 
certainly far less important than galaxy metallicity. 

These observations do not agree with the suggestion 
of Mannucci et al. (2011) "that the difference with the 
mass-metallicity relation is due to higher than average 
SFRs [of LGRB hosts] and that LGRBs with optical af- 
terglows do not preferentially select low-metallicity hosts 
among the star-forming galaxies." While the average 
SFR of LGRB hosts is indeed higher than that of typi- 
cal SDSS galaxies, this is because LGRBs do not choose 
galaxies based on number but rather based on SFR (as 
well as metallicity). The star- formation distribution of 
the LGRB hosts population tracks that of similar metal- 
licity SDSS galaxies, Indeed, due to the fact that the 
star-formation distribution of galaxies in the SDSS is 
largely independent of metallicity, they track the star- 
formation distribution of the entire SDSS as best as can 
be determined with only the 11 LGRBs for which we 
have good SFRs. The LGRB hosts population is explic- 
itly concentrated in the low metallicity end of the avail- 
able star-formation. While the LGRB hosts themselves 
may remain co nsistent with the mass , metallicity, and 
SFR relation of |Mannucci et aL (2010), this relationship 
is not sufficient to expiain the observed concentration of 
LGRBs in low metallicity star-formation. 

Nonetheless, the preference of LGRBs for low- 
metallicity hosts is not absolute. Three cases of LGRBs 
in roughly solar metallicity hosts (LGRBs 051022, 
020819B, & 050826) demonstrates that LGRBs can oc- 
cur at high metallicity. However, such events are quite 
rare. These high-metallicity objects are all at substan- 
tial redshifts (99% of the galaxies in the SDSS sample 
lie closer than the closest of these LGRBs), and thus the 
search volume needed to find them was large. As a result, 
our sample overemphasizes high-metallicity LGRBs com- 
pared to the distribution likely to be found in a volume 
limited survey. The presence of these objects does not 
substantially detract from our main conclusion, that on 
the whole, LGRBs prefer low metallicity. However, the 
absence among the three high metallicity LGRBs for any 
apparent preference for a lower metallicity environment, 
compared with the general star forming galaxy popula- 
tion of comparable brightness & redshift, may be a hint 
that the differential preference of LGRBs for lower metal- 
licity environments is reduced at high metallicities. We 
will discuss the implications of this in an up coming paper 
on high metallicity LGRBs ( |Graham|2012| 



m prep) 



It is generally accepted that the emission from LGRBs 
is relativistically beamed from collimated outflows and 
that rotation is required to produce the collimation 



dMeszaros fc Rees||1997||Sari||1999||Quilligan et al.||2002 
M cBreen et al.||2002l |jSlagataki et ;' al \\'M6\ |Toma et al'| 



Komissarov et al.|2009[[Bucciantini et al.|2009t|Liu et al. 
2010| |Nagakura et al.[|2011| [Yuan fc Zhang||2012| ). As 
increased metallici ty produces mass loss through stellar 
winds (c.f. | Nugis fc Lamers | |2000| |Crowther et al.||2002 " 



Vink fe deT^oter||2005| ), and these winds carry off an- 



gular momentum, it has been proposed ( Yoon fc Langer| 
2005 |Langer fc Norman |2006 ) that the observed preter- 
ence for low metallicity LGRB environments is a result 
of the need to ensure rapid progenitor rotation. However 
the predicted metallicity rang e of 1/10*' 1 or 1/100*' 1 solar 
required to limit mass loss does not match where obser- 
vations find a large drop in LGRB occurrence, at about 
a half to one-third solar metallicity. 

Any complete theory of LGRB formation must al so ac- 
count for the occasiona l LGRB at high metallicity. |Pod-| 
siadlowski et al. (2010) attempts to span this dichotomy 
through a binary common envelope ejection model in 
which low-metallicity helps drive mass transfer, but nu- 
clear burning during a merger phase provides the energy 
to eject the stellar atmosphere. However, this model also 
requires a much greater difference in metallicity (than 
present observationally) to yield an order of magnitude 
difference in the LGRB formation rate. This model's re- 
quirement of hydrogen and helium layer ejection does ex- 
plain the specific association of LGRBs with only Type Ic 
SNe but not the additional constraint that the SNe have 
the ultra fast ejecta of broad-line SNe events. Again, 
we defer a more detailed discussion of high metallicity 
LGRBs and their implications upon theory to our up- 
coming paper which specifica lly looks at high-metallicity 



LGRB hosts (Graham 



2012 in prep). 



2007| IKomissarov fc Ba rkov 2007| |Marscher et al.||2008 



We have shown that LGRBs have a strong intrinsic 
preference for low metallicity environments. Although 
nearly all LGRBs are thought to be associated with 
Type Ic-bl SNe, we have also shown that any preference 
for low-metallicities in the general Type Ic-bl popula- 
tion is minor in comparison. While LGRBs are clearly 
associated with star- format ion, they are highly biased 
by metallicity. Any use of LGRBs as tracers of star- 
formation must take this into account, though the fact 
that the production of LGRBs appears to turns on in 
force below one-half solar may give hope to those who 
wish to use LGRBs as tracers at high-redshift, where 
metallicities this low are the norm rather than the ex- 
ception. 

A number of people deserve great thanks for the help 
they have pr ovided us in doing this work. Lisa Kewley 
provided the |Kobulnicky fc Kewley ( 2004[) scale metal- 
licity code. Jarle Brinchmann provided Trcmonti ~et al.| 
(2004) scale metallicities for LGRB and SNe objects to 
match the SDSS object values in the MPA-JHU data 
products (though results in this metallicity diagnostic are 
not presented in this paper it provided a useful crosscheck 
throughout our analysis). Greg Aldering & Saul Perl- 
mutter provided unreleased proprietary SNFactory data 
for 5 Type II SNe. Kuntal Misra provided slit loss esti- 
mates for the LGRB population based on HST imaging 
data. Kuang-Han Huang provided multi-band SExtrac- 
tor photometry on Type Ic-bl SNe from SDSS data. In 
addition to these individuals, Rebekah Hounsell, Norbert 
Langer, Mario Livio, Selma de Mink, Maryam Modjaz, 
Colin Norman, and Rosie Wyse also provided many use- 
ful discussions. We thank them all. 
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APPENDIX 

A. METALLICITY GRADIENT FITTING AND APPLICATION 



LGRBs generally occur in the brightest regions of their host galaxies (Fruchter et al. 2006) and have small enough 
hosts that the majority of the total host light is observed via slit spectroscopy. SJNe however are often found in larger 
galaxies where there is a measurable metallicity difference between the location within the galaxy where the event 
occurred and th e galaxies center. Sin ce in most cases we lack observations of both metallicities, we follow the general 
methodology of |Modjaz et al. (2008) and employ an expected metallicity gradient to estimate the metallicity at SNe 
location when it was not specifically measured (see Figure [3]). In this paper, to better compare the physical properties 
of the host galaxies, we also reverse this methodology to estimate the central galaxy metallicity based on the SNe 
location metallicity and gradient when an independent R23 core galaxy metallicity measurement is not available (see 
F igure [8]). Here , we describe how we derive the galaxy metallicity gradients used for these estimates. 

iGarnett et al. ( 1997 1 Table 7 provides a compilation of metallicity gradients from a sample of thirty- two nearby spirals 
(see their Table 7). However, the provided metallicities are calculated via the Z94 ( Zaritsky et al.||l994[ ) methodology. 
W hile metallicity values ca n be converted from the Z94 scale to the KK04 scale used here via t he transform provided 



Kewley & Ellison (2008), this transform is for metallicity values not gradients. Fortunately Garnett et al. (19971 



also provides the central metallicites and the separation over which the gradient was observed. Thus we back compute 
the Z94 metallicites at the given distance, convert both these and the central metallicity values to the KK04 scale 
(via the Kewley fc Ellison]|200 8[ transform) and then recalculate the metallicity gradients in the KK04 scale. These 
gradient values are plotted as a function of Ms values in Figure [9] 
We then fit the metallicity gradient values (VZ) with the provided Mb values generating the function given below: 



VZ = -0.028 x M B - 0.60 ± 0.016 (dex/kpc) 



(Al) 



We note that beyond M B values brighter then approximately -21.5, this function extrapolates to an inverted gradient 
with galaxies increasing in metallicity with radius. To avoid this, we adopt a zero gradient for galaxies brighter than - 
21.5 (in practice this affects SN 1998ey only). We use this equation to estimate central metallicites from site metallicities 
and vice versa in the obvious manner. We estimate the metallicity gradient from the host Mb magnitudes, multiply 
by the separation distance to determine the metallicity shift and use this differential to convert between progenitor 
site and galaxy core metallicities as needed. R23 metallicites are used where available and in a cases where only an 
[N II] /Ha central metallicity value is available for the host and an R23 values for the SNe site the shifted R23 value is 
used for the central metallicity due to the respectively higher accuracy of the R23 diagnostic. 
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Figure 9. Metallicity gradient values vs. B band ab solute gala xy lumin osity. Gradient values calculated from the given me tallicity values 



-0.028 x M B - 0.60. 



performed 



yields Vmetallicity (dex/kpc) 



We note that this approach does not compensate for metallicity grad ient differences as a f unction of metallicity and 
attempted to correct for this by fitting the metallicity gradient values of Garnett et al.| ( |1997| as a function of both M b 
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and central metallicity (and using iteration when solving for the central metallicity) . However this proved problematic. 
We suspect the primary causes may be the small number of data points in the Garnett sampl e, the general di f ficulty 



of accurately transforming between metallicity indicators, and potentially problems with the Zaritsky et al. (19941 
itself, most particularly its failure to include an estimate of i onization. A more ac curate approach and an interesting 
exercise in itself would be to repeat the gradient analysis of Garnett et al. (19971 with the more modern metallicity 
diagnostics directly and then attempt a fit in Mg and metallicity however this exceeds the scope of this paper. 

Due to our inability to adjust our metallicity gradient calculations as a function of galaxy metallicity, we suspect 
that in cases with high SNe site metallicity (i.e. SN 1997ef) we are likely over estimating the central metallicity. It 
is apparent in the scatterplots (figures |3jJ3 and [8]) that while the luminosity metallicity relation is linear over most if 
its range for the brightest objects the sDSS population appears to flatten in metallicity at about ~2 solar. This is 
to be expected as the metallicity will eventually reach a maximum yield value, in terms of the metallicity gradient 
within the galaxies, as the central galaxy metallicites truncate at the yield while the enrichment continues in the 
remainder of the galaxy the metallicity gradient will progressively flatten from the center of the galaxies outward. 
This effect is progressively more noticeable for larger galaxies as the enrichment reaches saturation at progressively 
increasing radii. We also note that the [N II] /Ha diagnostic seems to allow metallicities beyond the maximum of those 
commonly returned by R23 of about log(0/H)+12 ^9.1, although the values below this remain comparable. This may 
be indicative of different yield values for Nitrogen vs. Oxygen or just an artifact of the different metallicity diagnostics. 

While an investigation of these effects is beyond the scope of this paper we do attempt a crude correction for these 
effects as follows. We limit the the measured [N II]/Ha metallicities to a maximum log(0/H)+12 value of 9.1, thus 
forcing it into agreement with the R23 range; however when downshifting the core metallicities to estimate a site value 
we begin with the actual [N II] /Ha value (in practice this is relevant to SN 1997dq only). We also place a variable 
limit on the maximum estimated upshifted central metallicities based on upon galaxy absolute magnitude. To do this 
we we determine the 90th percentile metallicity of the SDSS galaxies as a function of galaxy absolute magnitude and 
use this as a celling for the gradient upshifted values. Additionally in the case of SN 2006nx the [N II] /Ha ratio (used 
to break the degeneracy between the R23 upper and lower branches) is of sufficiently high error to not firmly exclude 
the lower branch solution thus the error bar on this object is extended to include this possibility. 
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Figure 10. Measured central galaxy (left) and event site (right) metallicity vs. B band absolute galaxy luminosity of LGRB (red squares) 
and Type Ic-bl SNe hosts (blue circles) without associated LGRB events (filled circles represent SNe selected in an untaragted manner 
whereas open circlers were from galaxy targeted SNe surveys and thus may be biased in galaxy properties by target selection) as well as 
star-forming galaxies from the SDSS (green small dots) and TKRS (teal diamonds) to provide a low and high redshift comparison sample 
respectively. Due to the absence of line flux error data on many objects a consistent err or of ±0.07 dex is assumed fo r R23 metallicities 
and ±0.15 for [N II] /Ha metallicities with all metallicities consistently recalculated in the IKobulnicky &: Kewleyl (|2004| scale as described 
in Section |3.1| Subsequently these plots are combined such that when both measured site and central galaxy values are not available the 
missing value is estimated from the other using the metallicity gradient methodology outlined in AppcndixlAko generate a complete set of 
site and central galaxy values as plotted in figures [3] & [8] respectively. LGRBs are assumed to occur in thebrightest regions of their host 
galaxies l |Fruchter et al.|2006[ l and used as both central and site values unless known otherwise (LGRBs 980425, 020819B, 060505). 
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B. APPARENT BIASES OF NON- TARGETED SNE POPULATIONS 

SNe found in surveys targeted at a population of galaxies are generally expected to be biased by the criteria used to 
select that galaxy population. In general, target surveys preferentially search larger & brighter galaxies. However, one 
might expect that SNe surveys that are not targeted but instead observe an area of sky without emphasis on particular 
galaxies would produce truly unbiased SNe populations. Here we find evidence that this is not the case. 

In Figure 11 we show that the untargeted surveys used by us appear to have a bias against finding SNe in high 
surface brightness hosts. We can obtain a crude estimate of the central surface brightness of the host galaxies by using 
the plug magnitudes from the SDSS survey. Given the SDSS fiber plugs are 3" in diameter, this corresponds to a 
physical scale of ~1.2 to 2.3 kpc across our volume limited sample. Thus for large galaxies which subtend the SDSS 
fiber, the fiber plug magnitudes can serve as a crude proxy for surface brightness. None of the non-targeted Type 
Ic-bl and only 2 out of 13 non-targeted Type II SNe have a g band plug magnitude brighter than 18. The latter ratio 
further reduces to 1 out of 10 non-targeted Type II SNe when only considering non-targeted SNe found in dedicated 
non-targeted SNe searches. Compared with the 48% of star-formation found in the volume-limited SDSS sample in 
galaxies with g band plug magnitudes brighter than 18 this suggests a strong, and we believe non physical, bias against 
SNe occurring in high surface brightness galaxies. This claim is only further strengthened by a detailed examination 
of the two exceptions. The higher surface brightness case (SN 2007es) is actually located on the outskirts of a large 
galaxy in a region with much less surface brightness than the core where the SDSS fiber observed. The other exception 
(SN 2001fb) was found not in a dedicated non-targeted SNe survey but by incidental o bservation presum ably involving 
an alert human rather than search software that often avoids the centers of galaxies ( Rest et al.||2001[ ). 

This surface brightness aversion could perhaps explain features seen in both our plots of cumulative population 
vs metallicity (Figure [5]) and SFR distribution (Figure p I. In the latter figure, the untargeted Type II SNe deviate 
strongly from the expected SFR distribution. (The small but predominately non-targeted set of Type Ic-bl SNe for 
which we have good SFRs also shows a similar but more muted bias). The non-targeted Type II discrepancy becomes 
quite alarming when we consider only events detected via dedicated non-targeted SNe searches, namely the SN factory 
& SDSS surveys. (Some non-targeted SNe are the result of incidental discoveries where the original target was a region 
of sky or a galaxy unrelated to the SNe host). That this bias in the non-targeted population is not even diminished 
when the analysis of Figure [6] is limited to similar mctallicities suggests that, at least for the Type II SNe, which is 
consistent with a surface brightness bias. This bias could also be involved in the tendency of the non-targeted samples 
to lie at slightly lower metallicities than the SDSS population. 

A bias against detecting non-targeted SNe in bight galaxies would also explain another curiosity; only a singe SNe 
in each of the non-targeted Type II and non-targeted Type Ic-bl SNe populations is in a host more luminous than 
-20 Mb (whereas about a third of the star- formation in our volume limited SDSS sample is in galaxies more luminous 
then -20 Mb)- Furthermore introducing a maximum luminosity cut at an Mb of -20 forces the non-targeted Type 
Ic-bl SNe and the general star- forming galaxy populations into good agreement. While our minimum luminosity cut 
at -18 Mb is justifiable given the completeness limit of the SDSS survey (at our volume limited redshift range), such 
a maximum luminosity would however be at best a crude compensation for our lack of plug magnitudes (or central 
surface brightness estimates) for the majority of our Type Ic sample. 
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Figure 11. SDSS g band fiber plug apparent magnitude vs. SFR for SNe. Since the SDSS survey uses a consistent fiber size the apparent 
fiber plug magnitude functions as a proxy measure for galaxy surface brightness. In the left plot the entire SDSS general star-forming 
galaxy sample meeting the criteria described within are plotted in the background. In the right plot the b ackground points are replaced 
with a one thousand object synthetic star- formation weighte d random sample a s described in Section |2.1.2| (The apparent duplication of 
two points in both the SDSS and Supernovae Factory with the Prieto et al. |2008 objects is exactly that - the same SNe was detected in both 
surveys). While the targeted sample (described in Section \2. 3. 2[ | shows good tracking of the synthetic p opulat ion the Supernovae Factory 
population (one of the two component samples in the non-targeted Type II SNe population - see Section |2.3.2| only shows this tracking in 
dimmer galaxies with only one of seven SNe host galaxies with a g band plug mag brighter then 18. A closer examination of this object 
shows a large galaxy where the SNe occurred on the outskirts in a region with much less surface brightness then the core where the SDSs 
fiber was located. (The SDSS SNe search population may also show this effect but with only three data points we cannot be certain). 
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